DISTRIBUTION  STATEMENT  A 

Approved  for  Public  Release 
Distribution  Unlimited 


i 


NATIONAL  ADVISORY  COMMITTEE 
FOR  AERONAUTICS  ' 

TECHNICAL  NOTE  2130 


CALCULATION  OF  TRANSONIC  FLOWS  PAST  THIN  AIRFOILS 
BY  AN  INTEGRAL  METHOD 
By  William  Perl 

Lewis  Flight  Propulsion  Laboratory 
Cleveland,  Ohio 

'  -t- 


Washington 
July  1950 


i 

'a 


Reproduced  From 
Best  Available  Copy 


Jtfflc  quality  msmmm 

fiq  r^o o-ib-sjyy 


20000801  149 


TABLE  OF  CONTENTS 


j  Page 

1 

SUMMARY . 1 

INTRODUCTION .  1 

COMPRESSIBLE  POTENTIAL  FLOW . .  •  • .  •  4 

Governing  equations  .....  .  4 

Small -perturbation  case .  8 

Prandtl  -G-lauert  range . 9 

Lower  transonic  range;  similarity  rules  .  10 

Explicit  solution  by  choice . of  curvature 

function, .  14 

Conditions  in  flow  field . 18 

Illustrative  examples  .  20 

Potential -limit  phenomenon  ...............  21 

LOWER  TRANSONIC  ASYMMETRIC  FLOW  WITH  TERMINAL  SHOCK  .  28 

Velocity  distributions . 28 

Terminal  shock  ......  .  33 

UPPER  TRANSONIC  FLOW  WITH  DETACHED  SHOCK . .  •  •  39 

Analysis . 39 

Some  comparisons  with  experiment .  46 

PRESSURE  DRAG;  USE  OF  MOMENTUM  INTEGRAL . 48 

CONCLUDING  REMARKS . 52  ' 

APPENDIX  A  -  SMALL -PERTURBATION  APPROXIMATION  OF 

CONTINUITY  INTEGRAL .  53 

APPENDIX  B  -  SMALL -PERTURBATION  APPROXIMATION  OF 

MOMENTUM  INTEGRAL .  57 

APPENDIX  C  -  TRANSONIC  SMALL -PERTURBATION  APPROXIMATION 

OF  OBLIQUE -SHOCK  RELATIONS . 60 

APPENDIX  D  -  TERMINAL -SHOCK  CONDITION  IN  LOWER 

TRANSONIC  RANGE.  . . 64 

APPENDIX  E  -  SYMBOLS . • . .  •  70 

REFERENCES  .  . . .  ♦  •  •  . .  72 


NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  2130 

CALCULATION  OF  TRANSONIC  FLOWS  PAST  THIN 
AIRFOILS  BY  AN  INTEGRAL  METHOD 
By  William  Pari 


SUMMARY 

A  method,  of  calculating  two-dimensional  compressible  flows  past 
thin  airfoils  is  presented  with  particular  reference  to  the  transonic 
speed  range.  The  method  is  based  on  the  integral  form  of  the  equa¬ 
tion  of  continuity  and  on  the  intrinsic  form,  in  terms  of  the  stream¬ 
line  curvature,  of  the  irrotationality  condition.  The  application  to 
the  transonic  range  conforms  with,  and  is  carried  out  within  the 
framework  of ,  the  transonic  similarity  theory.  The  results  include 
velocity  distributions  for  various  symmetric  sections  in  the 
continuous-potential-flow  speed  range  and  for  the  symmetric  biconvex 
airfoil  at  zero  angle  of  attack  through  the  transonic  speed  range 
and  the  associated  variation  of  terminal-  and  head-shock  location 
and  pressure  drag  with  Mach, number .  Comparisons  with  other  theories 
and  'with  available  experimental  data  indicate  at  least  qualitatively 
good  agreement.  The  location  and  the  movement  of  the  terminal  shock 
in  the  lower  transonic  speed  range  are  discussed,  as  well  as  the 
limiting-line  phenomenon. 


INTRODUCTION 

Considerable  effort  has  been  expended  in  recent  years  on  the 
investigation  of  the  flow  phenomena  that  occur  when  aerodynamic 
bodies  travel  at  transonic  speeds.  Experimentally,  the  main  fea¬ 
tures  of  steady  transonic  flow  patterns  about  airfoil  sections  have 
been  established.  (See,  for  example,  references  1  to  5.) 

Thus,  the  succession  of  flow  patterns,  of  which  an  approx¬ 
imate  theoretical  description  is  attempted  herein,  may  be  briefly 
outlined  as  follows :  Consider  the  zero-lift  flow  relative  to 
an  airfoil  that  is  symmetric  with  respect  to  its  chord  line 
(fig.  1).  At  low  subsonic  free-stream  velocities,  the  flow 
pattern  is  likewise  symmetric  with  respect  to  the  chord  line  and  is 
everywhere  subsonic.  Above  a  certain  subsonic  free-stream  velocity, 
a  symmetric- type  local  supersonic  region, bounded  by  the  locus  of 
sonic  speed  indicated  by  1  (fig.  1),  starts  to  develop  in  the 
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neighborhood  of  the  maximum- thickness  location  of  the  airfoil*  If 
the  airfoil  had  fore-and-aft  symmetry  with  respect  to  its  midchord 
location,  these  local  supersonic  regions  would  likewise  have  fore- 
and-aft  symmetry  with  respect  to  midchord.  As  the  free-stream 
velocity  is  further  increased,  the  fore-and-aft  type  of  symmetry 
of  the  flow  pattern  disappears  and  an  asymmetric  pattern  develops. 

The  local  supersonic  region,  indicated  by  the  sonic  boundary  2  in 
figure  1,  is  now  of  an  asymmetric  type  (that  is,  would  be  asymmetric 
fore  and  aft  even  if  the  airfoil  were  symmetric  fore  and  aft)  and 
is  terminated  on  the  downstream  side  by  an  approximately  normal 
shock.  Near  the  airfoil  especially,  this  terminal  shock  is  actually 
a  more  or  less  complex  shock  structure  depending  on  the  Reynolds 
number  and  interacts  with  the  boundary  layer  in  a  complex  manner. 

As  the  free-stream  velocity  increases  toward  sonic,  the  asymmetric- 
type  local  supersonic  region  increases  in  extent  and  the  terminal 
shock  increases  in  intensity  and  moves  downstream  toward  the  trail¬ 
ing  edge.  At  slightly  supersonic  free-stream  velocities,  an  addi¬ 
tional  shock,  the  bow  wave,  appears  ahead  of  the  airfoil  and 
approaches  it  as  the  free-stream  speed  is  increased.  The  subsonic 
part  of  the  flow  field  has  now  become  finite  in  extent,  as  indicated 
by  the  sonic  boundary  3  (fig.  1).  The  terminal  shock  has  reached, 
and  remains  in,  the  neighborhood  of  the  trailing  edge  and  begins  to 
slope  in  a  downstream  direction.  At  a  certain  supersonic  free-stream 
speed  the  forward  shock  attaches  at  the  nose  of  the  airfoil,  if  the 
nose  is  sharp-edged,  and  shortly  thereafter  the  flow  field  becomes, 
usually,  everywhere  supersonic,  except,  of  course,  in  the  boundary 
layer. 

On  the  theoretical  side,  much  of  the  effort  expended  on  the 
transonic  problem  has  gone  into  the  solution  of  the  mathematically 
well-defined  problem  of  the  continuous  compressible  potential 
flow  past  an  airfoil,  corresponding  to  the  symmetric- type  flow 
pattern  described.  (See,  for  example,  references  6  to  20.)  Emmons 
and  Maccoll  and  Codd  (references  18  to  20),  however,  have  included 
shocks  in  some  of  their  relaxation  calculations.  Progress  with  the 
methods  cited  has  so  far  not  been  rapid,  primarily  because  of  mathe¬ 
matical  and  computational  difficulties.  Perhaps  the  most  significant 
theoretical  result  that  has  emerged  is  the  probability  that  a  sub¬ 
sonic  free-stream  speed  exists  for  a  given  aerodynamic  shape,  such 
that  continuous  symmetric-type  flow  patterns  containing  local  super¬ 
sonic  regions  can,  at  least  in  principle,  be  derived  for  this  shape 
below  this  speed  but  not  above.  In  regard  to  the  other  main 
ingredient  of  real  transonic  flows,  namely,  the  boundary  layer  and 
its  interaction  with  the  terminal  shock,  interpretation  and  analysis 
is  only  in  a  preliminary  stage.  (See,  for  example,  references  3, 

4,  and  21  to  24.)  It  must  be  noted,  however,  that  progress  with  a 
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complete  theory  is  hampered  by  the  practical  difficulties  involved 
in  carrying  out  good  experiments  in  the  transonic  speed  range. 

In  the  present  report  an  integral  method  is  applied  to  the  cal¬ 
culation  of  transonic  flow  patterns.  Starting  from  the  equation  of 
continuity  in  integral  form  and  conservation  of  energy  in  the  inte¬ 
grated  form  of  Bernoulli's  equation,  which  suffice  for  the  analysis 
of  quasi  one- dimensional  flows,  two-dimensional  effects  are  then 
taken  into  account  by  inclusion  of  the  equation  of  motion  normal  to 
the  streamlines,  expressed  in  terms  of  the  streamline  curvature. 

By  an  appropriate  choice  of  the  streamline  curvature,  the  flow 
pattern  can  then  be  calculated  in  a  comparatively  simple  manner. 

This  method  has  long  been  known  and  applied  to  the  calculation 
of  flows  through  curved  channels,  for  example,  between  turbine 
blades,  or  in  wind  tunnels  (references  25  and  26).  Application  has 
recently  been  made  to  isolated  airfoils  in  continuous  compressible 
potential  flow  (references  27  to  29). 

The  small -perturbation  form  of  the  method  of  reference  29  is 
developed  herein.  The  resulting  simplifications  make  feasible  the 
approximate  calculation  of  some  of  the  more  complicated,  shock- 
containing,  flow  patterns  previously  described.  The  small -perturbation 
form  also  shows  explicitly  the  agreement  of  the  method  of  reference  29 
with  the  more  general  transonic  similarity  theory  of  von  Karman 
(reference  30)  and  Oswatitsch,  to  which  any  nonviscous-type  explicit 
solution  for  thin  airfoils  in  the  transonic  speed  range  must  conform. 

The  mathematical  simplicity  of  the  present  method  preserves  the 
possibility  of  interpreting  results  and  difficulties  at  all  stages 
in  terms  of  the.  basic  conservation  laws.  '  The  insight  thus  permitted 
into  the  nature  of  transonic  flows  makes  it  easy  to  understand,  for 
example,  the  previously  mentioned  phenomenon  of  an  upper-limiting 
subsonic  flight  speed  for  the  existence  of  continuous  potential 
flows . 

In  order  not  to  unduly  disturb  the  continuity  of  presenta¬ 
tion  of  the  main  ideas,  certain  mathematical  details  are  relegated 
to  appendices  A  to  D  and  a  list  of  the  main  symbols  used  is  con¬ 
tained  in  appendix  E. 

The  main  developments  of  this  report  were  worked  out  at 
Columbia  University  in  the  winter  of  1947-1948  under  the  stimulus 
of  a  lecture  given  by  Professor  Theordore  von  Karman  on  the  subject 
of  transonic  similarity.  The  report  was  completed  at  the  NACA  Lewis 
laboratory. 
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The  author  wishes  to  express  his  thanks  to  Professor  Theodore 
von  Karman  for  having  inspired  the  present  treatment  of  the  tran¬ 
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von  Karman  has  probably  done  most  to  demonstrate  the  power  and 
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the  present  theory  is  built,  namely,  the  integral-method  approach 
and  the  similarity  concept.  It  is  to  be  hoped  that  the  present 
structure  does  not  rest  too  lightly  on  the  foundations. 


COMPRESSIBLE  POTENTIAL  FLOW 

Governing  equations.  -  The  flow  is  considered  steady,  two- 
dimensional,  isentropic,  irrotational,  and  symmetric  with  respect 
to  the  airfoil  chord  (x-axis).  The  equations  of  motion  can  be 
expressed  in  terms  of  the  local  fluid  density  p,  speed  v,  pres¬ 
sure  p,  Mach  number  M,  the  corresponding  free-stream  values 
(denoted  by  the  subscript  o),  the  specific-heat  ratio  7,  and 
the  streamline  curvature  C,  as  follows  (fig.  2):  (The  main  symbols 
are  defined  in  appendix  E.) 

Equation  of  continuity: 


rv, 


nn 


p0vody 


pvdn 


(1) 


Irrotationality  condition  (reference  31) : 

s  +  Cv  -  0  <2> 

Bernoulli-state  equation: 


1  1 


The  differentials  dy  and  dn  in  equation  (l)  are  elements  of 
length  along  a  potential  line  in  the  free  stream  far  ahead  of  the 
airfoil  and  along  a  potential  line  in  the  region  of  the  airfoil, 
respectively.  The  lower  limits  of  integration  lie  on  the  zero 


1255 


NACA  TN  2130 


5 


streamline,  that  is,  on  the  x-axis  and  on  the  airfoil  contour.  The 
upper  limits  lie  on  a  streamline  DBF,  which,  as  yQ  approaches 

infinity,  becomes  parallel  to  the  x-axis.  It  is  assumed  for  the 
present  that  as  yQ  approaches  infinity  the  displacement  component 
6  in  the  y-direction  between  any  two  points  on  the  same  streamline, 
such  as  D  and  B,  approaches  zero.  This  assumption  is  valid  for 
circulation-free  (zero-lift)  flows  that  are  subsonic  at  infinity 
and  that  have  no  wakes  or  unclosed  portions  of  airfoil  contour 
extending  to  infinity  in  the  x- direction,  inasmuch  as  such  flows 
can  be  represented  at  infinity  by  a  collection  of  incompressible- 
type  source-sink  doublets  in  the  finite  part  of  the  plane,  and  for 
such  doublets  the  assumption  is  valid.  This  assumption  can  also  be 
regarded,  perhaps  more  generally,  as  restricting  the  flow  pattern 
to  that  produced  in  a  closed,  straight-walled,  two-dimensional  wind 
tunnel  in  the  limiting  case  in  which  the  walls  are  infinitely  for 
apart . 


The  element  of  length  dn  in  equation  (l)  is  now  approximated 
by  the  element  of  length  dy  in  the  y-direction,  that  is,  the 
potential  lines  are  replaced  by  the  lines  x  =  constant.  This 
replacement  is  exactly  valid  for  the  potential  line  starting  from 
the  midchord  location  of  a  symmetric  airfoil  that  has  also  fore- 
and-aft  symmetry.  In  appendix  A,  this  replacement  is  shown  to  be 
justified,  in  general,  for  the  case  of  small  perturbation  of  the  free- 
stream  flow  by  an  airfoil. 

Next  let  the  lower  limit  on  the  right  side  of  equation  (1)  start 
from  a  point  on  the  airfoil  at  which  the  airfoil  ordinate  is  Y. 

Set  p0,  vQ,  and  pQ  equal  to  unity,  so  that  the  corresponding 

local  quantities  p,  v,  and  p  henceforth  represent  nondimens ional 
fractions  of  free-stream  density,  velocity,  and  pressure,  respec¬ 
tively.  Furthermore,  let  the  unit  of  length  be  the  airfoil  chord  c 
so  that  the  lengths  x,  y,  and  l/c  henceforth  represent  non- 
dimensional  fractions  of  the  airfoil  chord  length.  Subtracting 
y0-Y  from  both  sides  of  equation  (1),  letting  yQ  approach  infin¬ 
ity,  and  using  the  preceding  assumption  for  6  and  approximation 
for  dn  yields 


Y  = 


(pv-l)dy 


Equations  (2)  and  (3)  become, respectively. 


,  dv 

d3r  ■  •  s 


(S) 
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The  form  of  the  continuity  equation  (4)  is  analogous  to  that  defining 
the  displacement  thickness  in  the  boundary- layer  theory.  The  airfoil 
ordinate  Y  here  takes  the  place  of  the  boundary- layer  displacement 
thickness.  Whereas  the  subtraction  and  limiting  process  leading  to 
equation  (4)  is  performed  in  the  boundary- layer  theory  mainly  for 
definiteness  in  defining  a  boundary- layer  thickness,  the  process  is 
here  essential  in  order  to  subtract  out  the  infinite  mass  flow 
crossing  a  potential  line  and  initially  introduced  perhaps  somewhat 
artificially  by  considering  the  fluid  as  moving  relative  to  the 
stationary  airfoil. 

The  streamline  curvature  C,  considered  positive  when  the 
streamline  curves  in  the  same  sense  as  the  corresponding  contour  of 
the  airfoil  at  the  maximum  thickness  location,  is  a  function  of  the 
position  coordinates  x  and  y  of  the  flow  field.  At  a  given 
chordwise  location  x,  the  streamline  curvature  is  a  function  of  y 
that  varies  continuously  from  the  curvature  Ca  of  the  airfoil 
contour  at  y  =  Y  to  the  value  zero  at  y  =  oo.  it  proves  conven¬ 
ient  to  consider  the  curvature  C  at  a  given  x  as  a  function  of 
the  local  velocity  v  rather  than  of  y.  In  the  applications  made 
herein,  v  will  be  a  monotonic  function  of  y.  Hence,  at  least  for 
such  applications,  C  can  be  expressed  as 

C  =  Caf(x,v)  (7) 

The  assumption  that  the  streamline  curvature  can  be  as  well 
expressed  by  a  function  of  x  and  v  as  by  a  function  of  x  and 
y  may  be  more  generally  justified  as  follows:  Of  the  four  basic 
variables  of  the  problem,  namely,  the  vector  velocity  given  by  the 
two  variables  v  and  0,  where  0  is  the  angle  of  inclination  of 
the  local  stream  direction  to  the  x-axis,  and  the  vector  distance 
given  by  the  two  variables  x  and  y,  two  variables  are  to  be 
regarded  as  independent  and  two  as  dependent.  The  solution  of  the 
so-called  direct  problem  of  airfoil  theory  is  required  in  the  form 
where .  v(x,y),  0(x,y),  and  in  particular  the  velocity  V(x,Y)  on 

the  given  airfoil  contour  0a(x,Y)  are  finally  known.  The  analysis 
of  the  problem,  however,  can  be  made,  at  least  in  principle,  with 
any  two  of  the  variables  considered  as  independent  and  the  remaining 
two  as  dependent.  The  hodograph  method,  for  example,  is  based  on  v 
and  0  as  the  independent  and  x  and  y  as  the  dependent  variables. 
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The  present  method  is  based  on  v  and  x  as  the  independent  and  0 
and  y  as  the  dependent  variables.  Conceivably,  other  possible 
combinations  might  be  advantageous  in  particular  problems.  Similar 
remarks  apply  with  the  potential,  or  the  stream  function,  or  both, 
replacing  one  or  two  of  the  present  variables. 

The  present  method  is  thus  a  sort  of  semihodograph  method. 
Whereas  it  retains  the  principal  advantage  of  the  hodograph  method, 
namely,  that  the  velocity  v  and  hence,  by  equation  (6),  the  density 
p  becomes  essentially  an  independent  variable,  it  avoids  the  well- 
known  disadvantage  of  the  hodograph  method;  namely,  difficulty  in 
satisfying  the  boundary  conditions  in  obtaining  flows  past  given 
fixed  airfoils  as  a  function  of  Mach  number. 

Equation  (4)  becomes,  upon  substitution  of  equations  (5)  and 

(7), 


in  which  Y  is  the  sought-for  velocity  at  the  airfoil  at  chordwise 
location  x,  the  quantity  x  being  held  constant  in  the  integration. 
Equations  (6)  and  (8)  constitute  an  expression  for  the  velocity  V 
at  a  point  on  an  airfoil  in  terms  of  the  shape  parameter  YCa  of  the 
airfoil  at  that  point,  the  free-stream  Mach  number  M0,  and  the 

curvature  function  f(x,v).  (From  the  point  of  view  of  the  inverse 
problem,  that  is,  designing  an  airfoil  to  satisfy  a  prescribed 
velocity  distribution,  the  shape  parameter  YCa  would  be  regarded 

as  dependent  upon  V,  Mg,  and  f.)  The  basic  reason  for  expressing 
the  airfoil  velocity  V  in  terms  of  a  streamline-curvature  function 
f  is  that,  in  the  form  of  equation  (8),  Y  is  relatively  insensitive 
to  choice  of  f .  As  will  be  seen,  a  reasonable  choice  of  f  can 
therefore  be  made  that  leads  with  comparatively  little  computation  to 
fairly  accurate  values  of  V. 


The  relation  between  the  local  velocity  v  and  the  lateral 
distance  y  at  a  given  x  is  obtained  by  integration  of  equation  (5) 
from  an  airfoil  point  (x,Y)  or  (x,V)  to  a  field  point  (x,y)  or 
(x,v)  as 


v 
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Solution  of  equation  (8)  for  V  enables  the  local  velocity  v(x,y) 
at  any  point  in  the  flow  field  to  be  determined  from  equation  (9). 

The  streamline  inclination  0(x,y)  can  then  also  be  obtained,  as 
will  be  subsequently  shown. 

The  only  approximation  thus  far  made  is  to  replace  the  element 
of  length  dn  along  a  potential  line  by  the  element  of  length  dy 
perpendicular  to  the  chord  line.  The  resulting  equations  (6),  (8), 
and  (9),  together  with  an  assumed  form  for  the  curvature  function 
f(x,v),  have  been  applied  in  reference  29  to  calculate  (finite 
perturbation)  compressible  potential  flows  past  airfoils.  An 
approximation  for  the  density  p(v)  is  next  made,  which  leads  to 
the  general  small- perturbation  form  of  the  preceding  theory. 

Small-perturbation  case.  -  By  small  perturbation  is  herein 
meant  that  the  local  velocity  increment  v-1  produced  by  the  airfoil 


is  small  relative  to  unity.  The  density  p  is,  by  equation  (6), 
thereby  expressible  as  a  power  series  in  v-1.  Thus,  the  numerator 
of  the  integrand  in  equation  (8)  can  be  developed  as 

p  -  I  =  pv  -  T  v2  +  .  .  . 

V  C 

(10) 

where 

v  sv  -  1 

(11a) 

5  1  -  Mq2 

(lib) 

_  M02  4 

Tcsl  +  -r-  2  M° 

(11c) 

urn  rc s  r=2£ 

(lid) 

Mo-»l 


The  velocity  increment  A  at  the  airfoil,  a  normalized  local- 
velocity  variable  z,  and  a  Mach  number  parameter  K  are  next 
defined  as 
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A  =  v  -  l 

„  _  V  v  -  1 

A  "  v  -  l 


is 


iL 


1-M 


(tD^  (tD 


O 

Tfl 


(12a) 

(12b) 

(12c) 


In  terms  of  the  reduced  local  velocity  z, 
(equation  (7))  can  be  expressed  as 

C  =  Ca  f(x,z,K) 


the  curvature  function 

(13) 


Equation  (13)  can  easily  be  shown,  from  the  general  transonic  sim¬ 
ilarity  theory  of  reference  30,  to  be  the  most  general  possible 
form  for  the  streamline- curvature  function  in  the  small-perturbation 
transonic  range.  (This  possibility  was  pointed  out  by  Professor 
Theodore  von  Kerman.)  Conversely,  it  will  be  presently  verified 
that  the  use  of  equation  (13)  with  the  present  formulation  of  the 
equations  of  motion  leads  to  the  transonic  similarity  laws  of  refer¬ 
ence  30.  The  curvature  function  f(x  z,K)  in  equation  (13) 
must  satisfy  the  boundary  conditions  f(x,l,K)  =  1  at  the  airfoil 
and  f(x,0,K)  =0  at  infinity. 


Substitution  of  equation  (10)  into  equation  (8)  and  changing 
the  variable  of  integration  from  v  to  z  by  equations  (lla), 
(12),  and  (13)  yields 


EcA  dz 

M  )  f(x,z,K) 


(14) 


Similarly,  the  relation  (9)  between  local  velocity  v  and  lateral 
distance  j  becomes,  for  small  perturbations, 


y 


J 


z 


dz 

f(x,z,K) 


(15) 


Prandtl-Glauert  range.  -  In  the  subsonic  or  Prandtl-Glauert 
range,  defined  by  0  <  Mq  <  <  1  or  K->  o>,  the  dependence  of  the 
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curvature  function  (equation  (13))  on  K  must  disappear.  Hence, 
in  the  small-perturbation  case,  A  <  <  1,  and  for  sufficiently 
subsonic  free- stream  Mach  numbers  MQ,  the  second  term  in  the 
integrand  of  equation  (14)  can  be  neglected  relative  to  the  first, 
and  equation  (14)  can  be  written  as 

YC 

-f  =  F(x)  (16) 

pA2 


where  F(x)  is  a  function  of  x  of  the  order  of  magnitude  of 
unity.  For  smoothly  curved,  thin  airfoils,  the  shape  parameter  YCa 
can  be  expressed  as  the  product  of  the  square  of  the  thickness  ratio 
T  (T  =  maximum  thickness/chord  length)  and  a  function  of  x  of  the 
order  of  magnitude  of  unity  over  most  of  the  airfoil.  Equation  (16) 
can  therefore  be  written  as 

A  =  -L  G(x)  (17) 

where  T  <<  1  and  G(x)  ~  1.  Equation  (17)  expresses  the  basic 
result  of  thin-airfoil  theory  that  the  velocity  increment  at  chord- 
wise  location  x  of  one  of  a  series  of  thin  airfoils,  all  having 
the  same  geometric  shape  but  differing  in  thickness  ratio  r,  is 
proportional  to  the  thickness  ratio  t  of  that  airfoil.  Equa¬ 
tion  (17)  also  expresses  the  Prandtl-Glauert  rule  that  the  effect 
of  compressibility  in  the  lower  subsonic  range  is  to  increase  the 

airfoil  velocity  increment  by  the  factor  1  /VI  -M02. 

Lower  transonic  range;  similarity  rules.  -  As  the  free- stream 
Mach  number  Mq  is  increased  toward  unity,  or  [x— >0,  the  airfoil 
contour  being  such  as  to  maintain  a  small  velocity  perturbation 
A  <  <  1,  a  value  of  \i  <  <  1  must  eventually  be  reached  where 
TA/\i  becomes  comparable  to  unity.  The  second  term  in  the  integrand 
of  equation  (14)  is  then  comparable  with  the  first.  Additional 
terms  in  the  expansion  of  the  integrand,  however,  evidently  would 
remain  negligible  relative  to  the  first  two  terms.  Hence,  at  tran¬ 
sonic  speeds  defined  by  ~  0,  instead  of  the  relation  of  equa¬ 
tion  (16)  between  the  two  parameters  x  and  YCa/pA^,  equa¬ 
tion  (14)  constitutes  a  relation  among  four  parameters;  thus 

k’  i2’  ^  m2^)"  ° 


H 


(18) 
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The  last  three  parameters  indicated  in  equation  (18)  are  not  inde¬ 
pendent,  however,  but  can  be  expressed  in  terms  of  two  independent 
parameters.  Thus,  using  the  third  parameter  of  equation  (18)  to 
eliminate  A  from  the  second  yields  the  equivalent  expression  of 
equation  (14) 


I 


■U 

(Tr)2/3'  v 


(19) 


or  also,  using  the  second  parameter  of  equation  (19)  to  alter  the 
third, 


J 


u  .  r1/3A\  _ 

FT75'  ~^)  " 


(20) 


Inasmuch  as  the  pressure  coefficient  P  at  a  point  on  the  airfoil 
is  obtained  in  terms  of  the  velocity  from  equation  (3)  in  the  small- 
perturbation  case  as  *  . 


-2A 


(21) 


relation  (20)  can  be  written  as 

t2/3 

P  -  pyj  9>  (x,K)  (22) 

The  function  in  equation  (22)  is  determined  by  the  shape  of  the 
airfoil  and  is  of  the  order  of  magnitude  of  unity,  Equations  (22) 
and  (12c)  agree  with  the  transonic-similarity  theory  as  regards 
velocity  distributions  on  airfoils. 

By  multiplying  equation  (15)  by  */“jT  and  using  the  preceding 
relations ,  there  results  for  the  velocity  increment  in  the  flow 
field 

(tiFVy),  k] 


V 


(23) 
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where  JS ^  and  J^2  are  functions  of  the  order  of  magnitude  of  unity 

and  W(y-Y) ,  the  affinely  transformed  lateral  distance  of  the 
Karman  theory  (reference  30)  is  of  the  order  of  magnitude  of  unity 
by  equation  (15).  Equation  (23)  agrees  with  the  transonic-similarity 
theory  as  regards  conditions  in  the  flow  field.  Equation  (23) 
reduces  to  equation  (22)  on  the  airfoil. 


Equation  (22)  states  that  if  the  pressure  coefficient  has  a 
certain  value  P^  at  a  chordwise  location  z  on  an  airfoil  of 

thickness  ratio  in  a  flow  at  speed  MQ  ^  ~  1  and  in  a  fluid 

of  specific-heat  ratio  y then  the  pressure  coefficient  P2  at 

the  same  chordwise  location  of  a  geometrically  similar  airfoil  of 
thickness  ratio  Tg  in  a  fluid  of  specific-heat  ratio  y  ^  and  at 

a  speed  MQ^2  given  by 

(24) 

(y2)2/3 

will  be  related  to  P-^  by 


p2 

Pi 


(25) 


The  Prandtl-Glauert  rule  is  included  in  equation  (22)  if  an  inverse 
square  root  is  chosen  for  the  K-dependent  part  of  the  function  . 
This  result  necessarily  follows  because  the  Prandtl-Glauert  case  is 
obtainable  from  the  general  formulation,  equations  (14)  and  (19), 
by  letting  rA/u->  0. 

A  further  similarity  rule  that  concerns  the  local  Mach  number 
may  be  noted.  In  isentropic  flow,  the  local  Mach  number  M  is 
related  to  the  local  velocity  v  and  the  free- stream  Mach  number 
MQ  as  follows : 

-  mV 

M2  =  - - -  (26) 

.  1  -  Zgi 

In  the  small-perturbation  case  v-l=  v  <<  1,  equation  (26) 
becomes 
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M^-l  =  -  |a  +  2rMV  (27) 

where 

rM  *  Uo  (>■  +  *o)  (28»> 

and 

iimrM=r  =  z+i  (28b) 

M0-*  1 

Equations  (23) ,  (27),  and  (28b)  yield  the  transonic  similarity  rule 
for  local  Mach  number 


■  —  ~  )-  (TH^fy-T),  K~|  (29) 

m2/ 

Equation  (29)  may  he  applied  in  a  manner  similar  to  equations  (24) 
and  (25)  to  relate  local  Mach  numbers  and  Mg  at  the 

corresponding  points 


X1  ~  x2 

<Tiri>1/3<vV  ■  <t/2>1/3<vV 


(30a) 

(30b) 


in  the  flow  fields  of  similar  airfoils  and  Tg  in  fluids 

characterized  by  T and  Tg  at  corresponding  free- stream  Mach 
numbers  and  MQ^g  given  by  equation  (24).  Instead  of  equa¬ 

tion  (25),  the  following  relation,  obtained  from  equation  (29), 
applies: 

Mg2-1  Mx2-1 

^^-■^73  <31> 

Various  special  cases  of  this  similarity  rule  have  been  previously 
given  in  references  4,  32,  and  33. 
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It  may  be  noted  that  although  the  use  of  T -  (y+l)/2  In  such 
expressions  as  equations  (23)  and  (29)  is  valid  in  the  transonic 
limit  Mq  =  1,  in  applications  for  vhich  MQ  is  finitely  different 

from  unity  equations  (28a)  and  (lie)  for  and  rc,  respectively, 

should  more  properly  be  used.  The  reason  is  that  in  the  approxi¬ 
mating  equations  (10)  and  (27),  only  the  small-perturbation  approxi¬ 
mation  V<<  1  was  made  and  not  the  transonic  approximation 
\i  <<  1.  This  transonic  approximation  was  made  in  equation  (18), 
and  so  forth,  only  to  show  the  agreement  between  the  common  tran¬ 
sonic  limiting  value  of  and  Tc  with  the  T  of  reference  30. 

Explicit  solution  by  choice  of  curvature  function,  -  Before 
choosing  a  form  for  the  curvature  function  f  =~c/ca  and  thence 

calculating  the  flow  pattern  explicitly,  a  few  general  remarks, 
which  may  be  helpful  in  judging  the  degree  of  approximation  of  the 
resulting  solution,  will  be  made*  The  present  formulation  of  the 
airfoil  flow  problem  is  in  terms  of  longitudinal  distance  x  and 
local  velocity  v  as  independent  variables  and  local  streamline 
curvature  C  and  lateral  distance  y  as  dependent  variables, 
that  is,  the  solution  for  given  stream  Mach  number  MQ  is 
required  in  the  form  C(x,  v),  y(x,  v).  Various  possible  pro¬ 
cedures  for  obtaining  this  solution  suggest  themselves.  For  example, 
an  iterative  procedure  would  be  to  assume  in  first  approximation  a 
curvature  function  C(x,  v),  and  thence  calculate  y(x,  v)  and, 
in  particular,  the  streamlines,  using  equations  (l)  to  (9).  The 
curvature  of  these  streamlines  would  then  yield  a  second  approxi¬ 
mation  C(x,  v)  and  the  procedure  would  be  repeated  until  satis¬ 
factory  convergence  was  achieved.  Again,  a  network-type  procedure 
would  be  to  set  up  the  continuity  integral  (l)  between  n  pairs  of 
streamlines  rather  than  between  the  one  pair  (airfoil  and  infinity) 
in  equation  (4).  The  number  n  would  increase  in  successive  steps 
of  the  calculation.  In  each  step  a  curvature  function  based  on  the 
curvature  values  obtained  in  the  preceding  step  would  be  used  to 
yield  by  simultaneous  solution  of  the  continuity  and  distance  equa¬ 
tions  (l)  and  (9)  an  increased  number  of  curvature  and  corresponding 
distance  values  for  use  in  the  next  step. 

The  rapidity  of  convergence  of  the  preceding  types  of  procedure 
could  be  expected  to  depend  strongly  on  the  initial  choice  of  cur¬ 
vature  function  in  the  first  approximation.  It  is  therefore  sensible 
for  computational  reasons  to  choose  an  initial  curvature  function  in 
the  light  of  the  information  supplied  by  experiment  and  other  theory 
as  to  the  behavior  of  the  streamline  curvature  in  the  type  of  flow 
field  under  consideration. 
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Only  the  first  approximation  is  considered  herein  and  free  use 
is  made  of  experiment  and  other  theory  in  setting  up  an  initial 
curvature  function.  In  view  of  the  preceding  remarks,  this  use  of 
experiment  and  other  theory  should  not  be  construed  as  inherently 
limiting  the  degree  of  rigor  in  principle  attainable  in  solutions 
by  the  present  approach.  This  ultimate  degree  of  rigor  can  be 
expected  to  be  about  the  same  as  with  existing  perturbation  or  net¬ 
work  formulations  of  the  airfoil  flow  problem. 

An  explicit  form  for  the  curvature  function  C/Ca  =  f(x,  z,  K) , 
equation  (13),  in  the  small -perturbation  case  is  now  set  up.  This 
function  must  satisfy  the  boundary  conditions  of  the  problem;  namely, 

at  airfoil, 

y  =  T 

v  =  A 

z  =  l 

c  =  ca 

f(x,l,K)  =  1 


)  (32a) 


at  infinity, 


y  — » eo 
V-40 
z  — »0 
C-*0 

f  (x,0,K)  —>  0 


'N 


> 


J 


(32b) 


Furthermore,  it  may  be  expected,  at  least  in  subsonic  flow  and  at 
the  chordwise  locations  near  the  maximum  velocity  or  thickness 
regions  of  the  airfoil  (which  are  of  most  interest),  that  the  stream¬ 
line  curvature  will  decrease  monotonically  from  the  positive  value 
Ca  at  the  airfoil  to  zero  at  infinity.  A  simple  curvature  function 
that  satisfies  these  requirements  is  (reference  29) 

f  =  zr  (33) 

in  which  the  parameter  r  is  a  function  of  x  and  K  to  be 
determined. 
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Substitution  of  equation  (33)  into  equation  (14)  and  integra¬ 
tion  yields 


YCa 


2 


ma A. 

2-r 


rcA3 

3-r 


(34) 


Equation  (34)  gives  the  local  velocity  increment  A  on  a  symmet¬ 
ric  airfoil  in  compressible  potential  flow  in  terms  of  the  local 
shape  parameter  YCa,  the  Mach  number  parameter  p,  the  specific- 
heat-ratio  parameter  rc,  and  the  parameter  r. 

In  order  to  determine  r,  it  is  first  noted  that  to  satisfy 
the  boundary  condition  at  infinity  (equations  (32b)),  r  must  be 
positive  and  the  integral  for  the  lateral  distance  y  in  equa¬ 
tion  (15)  must  diverge  at  z  =  0,  and  therefore  r  >1.  In 
addition,  the  continuity  integral  (equation  (14))  must  converge  at 
z  =  0,  which  it  does  for  r  <  2.  Hence  r  is  restricted  to  the 
range 


1  <  r  <  2  (35) 

In  the  speed  range  for  continuous  potential  flow,  the  flow  pattern 
is  largely  subsonic,  as  mentioned  in  the  INTRODUCTION.  Hence,  in 
this  speed  range  the  determination  of  r  will  be  completed  (l)  by 
assuming  r  to  be  independent  of  K,  and  (2)  by  requiring  equa¬ 
tion  (34)  to  yield  the  (presumed)  known  velocity  distribution 
Ai(x)  in  the  incompressible  case  Mc  =  0.  Conditions  (l)  and  (2) 

are  appropriate  to  the  Prandtl-Glauert  speed  range.  In  the  higher 
speed  ranges  discussed  in  subsequent  sections  in  which  the  flow 
pattern  differs  very  substantially  from  Prandtl-Glauert  type,  these 
conditions  on  r  will  be  discarded  for  new  ones.  By  condition  (2) 
and  equations  (lib)  and  (lie),  equation  (34)  yields  for  the  deter¬ 
mination  of  r 


YCa 


2-r  3-r 


(36) 


The  values  of  r  determined  by  equation  (36)  turn  out  to  be  in  the 
range  of  equation  (35)  over  the  maximum- velocity  region  of  the 
airfoil,  for  which  the  curvature  function  (equation  (33))  was  set 
up.  This  type  of  inner  consistency  will  recur  often  in  this  report. 

The  second  term  on  the  right  side  of  equation  (36)  can,  in 
principle,  be  neglected  in  the  small-perturbation  case  considered. 
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Better  results  are  obtained,  however,  in  applications  to  airfoils 
of  finite  thickness  ratio  (of  the  order  of  0.1)  if  it  is  retained. 
Some  theoretical  explanation  for  this  improvement  and  hence  retention 
of  the  term  in  question  can  be  made  as  follows:  The  right  side  of 
equation  (34)  actually  consists  of  the  first  two  terms  of  an  infinite 
series  in  the  velocity  increment  A,  which  is  evident  for  any  cur¬ 
vature  function  f,  from  equations  (8),  (10),  and  (14).  If  the 
coefficients  multiplying  the  powers  of  A  are  of  comparable  mag¬ 
nitude,  as  in  equation  (36),  it  may  therefore  be  expected  that  for 
small  but  finite  A  the  first  two  terms  of  the  series  would  be  a 
better  approximation  to  the  complete  series  than  the  first  term 
alone. 


Thus,  finally,  in  order  to  calculate  a  velocity  increment  A 
at  a  chordwise  location  x  on  a  symmetric  airfoil  at  zero  lift  at 
Mach  number  MQ  by  the  present  method,  the  shape  parameter  YCa 

and  the  low- speed  velocity  increment  Ai  at  that  point  are  deter¬ 
mined  first.  The  quadratic  equation  (36)  for  r  is  then  solved 
for  the  value  of  r  in  the  range  of  equation  (35).  (if  such  a  root 
of  equation  (36)  does  not  exist,  the  present  explicit  solution  based 
on  the  curvature  function  (equation  (33))  fails.)  While  keeping 
this  value  of  r  constant  as  Mq  varies  and  using  the  values  of  p 

and  Tc  corresponding  to  the  desired  Mq,  equation  (34)  is  solved 

for  the  desired  velocity  increment  A(x)  at  free-stream  Mach  number 
Mc.  For  convenience  in  solving  the  cubic  equation  (34)  for  A  for 

various  values  of  the  other  parameters,  these  other  parameters  can  be 
combined  as  follows:  Let 


g/'2~r\  A 

x  a  12^  tX 

(3-r)2  p3 

Equation  (34)  then  becomes 


(37a) 

(37b) 


(38) 


Equation  (38)  is  plotted  in  figure  3.  The  most  significant  feature, 
the  double- valuedness  of  (positive)  A  for  given  X,  vill  be  sub¬ 
sequently  discussed  in  connection  with  the  potential- limit 
phenomenon. 
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Conditions  in  flow  field,  -  Once  the  velocity  increment  A  at 
the  airfoil  has  heen  obtained,  the  variation  of  local  velocity  v 
in  the  field  vith  distance  y  from  the  airfoil  can  be  determined 
from  equation  (15).  With  the  curvature  function  (equation  (33)), 
equation  (15)  yields 


-1 

(r-l)(y-Y)Ca]r"1 

A 


(39) 


The  streamline  curvature  is  expressed  in  terms  of  y  by  equa¬ 
tions  (33)  and  (39)  as 


f 


(r-1) (y-Y)C 


1  + 


-r 

r-1 


(40) 


The  preceding  equations  show  that  both  the  local  velocity  and  the 
streamline  curvature  decrease  monotonically  vith  y  from  their 
values  at  the  airfoil  to  zero  at  infinity.  The  actual  manner  of 
this  decrease  depends ,  at  a  given  chordwise  location  x,  on  the 
airfoil  curvature  Ca  and  on  the  equations  of  motion  through  the 

parameters  r  and  A.  Actually  there  is  therefore  not  a  great  deal 
of  arbitrariness  in  the  values  of  streamline  curvature  obtained  here* 
Once  a  basic  form  such  as  equation  (33)  is  chosen,  the  actual  values 
of  streamline  curvature  in  the  flow  field  are  determined  by  the 
equations  of  motion  and  in  such  a  way  as  to  satisfy  these  equations 
on  the  average  or  in  the  large  (inasmuch  as  the  integral,  rather 
than  the  differential,  form  of  the  continuity  equation  was  used). 

The  insensitivity  of  airfoil  velocity  and  hence  also  of  the  flow 
pattern  generally  to  the  choice  of  a  particular  form  of  curvature 
function  such  as  equation  (33)  is  thus  plausible. 

A  streamline  in  the  flow  field  can  be  traced  out  as  follows: 

A  streamline  is  defined  by  a  constant  value  of  yQ  in  the  continuity 

equation  (l),  which  can  be  written  as 


y 


o 


(41) 
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Eliminating  v  in  equation  (41)  by  the  irrotationality  relation 
(5)  gives  the  equation  of  the  streamline 

y  =  S(x,y0)  (42) 

Thus  jQ  is  in  effect  the  stream  function.  In  the  small- 
perturbation  case,  equation  (4l)  can  be  transformed  as  follows: 


Elimination  of  z  in  equation  (43)  by  use  of  equation  (39)  gives 
the  equation  of  the  streamline  (equation  42).  The  inclination  6 
of  the  local  stream  direction  to  the  x-axis  can  then  be  obtained 
as  the  slope  of  the  streamline. 


It  is  of  interest  to  note  that  for  r  =  3/2,  which  is  in  the 
middle  of  its  permissible  range  (equation  (35)),  the  local  velocity 
increment  far  from  the  airfoil  is  given  in  the  Prandtl-Glauert 
range  by  equations  (39)  and  (16)  as 


lim  v 
y->  » 


yV^ 


(44) 


This  expression  agrees  with  that  obtained  for  a  source- sink  doublet 
by  conventional  methods. 
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Illustrative  examples,  -  In  order  to  illustrate  the  preceding 
explicit  solution,  consider  first  the  Kaplan  section  (reference  6) 
of  thickness  ratio  T  =  0.1.  This  contour,  indicated  in  figure  4, 
can  he  expressed  in  parametric  form  as  (reference  34) 


2Y  = 


n  =  I  t(s1 


sincp  -  —  sin  3cp 

\  3 


(45a) 


2x-l  a  l 


cos  cp  +  ± 


cos  3cp 


(45b) 


where  0  <9<2n.  When  the  shape  parameter  YCa(x),  as  deter¬ 
mined  by  equation  (45),  and  the  incompressible  velocity  distribu¬ 
tion  A^(x),  as  determined  by  conformal  transformation,  are  sub¬ 
stituted  into  equation  (36),  the  parameter  r(x)  is  found  to  vary 
from  r  =  1.301  at  the  midpoint  |  =  0  to  r  =  2  at  (=  0.645 
(because  =  0  but  YCa  ^  0  at  £  =  0.645).  Thus,  by 
equation  (35),  the  preceding  explicit  solution  fails  outside  the 
range  0  <  |  ||<  0.645.  The  velocity  distributions  A(x)  at  var¬ 
ious  Mach  numbers  Mg,  as  obtained  from  equation  (34)  within  this 

chordwise  range,  are  shown  in  figure  4.  For  comparison,  the  results 
obtained  (reference  29)  without  making  the  small-perturbation  approx¬ 
imation  of  equation  (10 )  for  the  density  and  the  results  of  Kaplan 
(reference  6)  are  shown.  The  results  of  Kaplan  were  obtained  by 
substituting  a  series  for  the  potential  in  powers  of  T  into  the 
partial  differential  equation  for  the  potential  and  successively 
solving  the  Poisson  equations  obtained  by  equating  orders  of  mag¬ 
nitude  in  T.  The  calculation  was  carried  out  to  the  third  order, 

3 

T  . 


The  three  sets  of  distributions  in  figure  4  are  seen  not  to 
differ  greatly  except  for  Mach  numbers  near  the  potential  limit, 
that  is,  the  maximum  Mach  number  for  which  the  method  yields  a 
solution.  This  limiting  case  will  be  discussed  in  the  next  section. 

As  a  further  illustration,  A  (Mg)  at  the  maximum  thickness 
locations  of  the  three  T=  0.1  sections  -  Kaplan,  symmetric 
biconvex,  and  ellipse  -  as  calculated  by  the  small-perturbation 
integral  method,  the  finite-perturbation  integral  method  (refer¬ 
ence  29),  and  by  various  other  authors  is  shown  in  figure  5. 

The  results  of  Hantzsche  (reference  7)  for  the  ellipse  were 
obtained  by  a  method  similar  to  that  used  by  Kaplan.  The  calculations 
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were  carried  to  the  third  order  in  a  thickness-ratio  parameter,  which 
in  successive  orders  equalled  T,  T2,  and  T"5  +  0(t^  log  t), 
respectively. 

In  the  method  of  Gothert-  and  Kawalki  (reference  17),  the  influ¬ 
ence  of  compress ibility  was  represented  by  a  finite  number  (12  per 
quadrant)  of  sources  and  sinks  (Rayleigh-  Janzen  idea)  located  in  the 
incompressible-flow  field,  as  corrected  by  the  Prandtl-Glauert  rule, 
outside  a  given  airfoil.  The  source-sink  strength  used  for  calcu¬ 
lating  velocity  increments  caused  by  compressibility  was  computed 
by  iteration  from  the  flow-field  velocities  themselves  so  as  to 
satisfy  continuity.  As  many  successive  steps  in  the  calculations 
( ~  10  near  the  potential  limit)  were  made  as  required  for  obtaining 
a  self-consistent  compressible-flow  field  that  could  be  regarded  as 
an  incompressible-flow  field,  with  the  given  spatial  arrangement  of 
sources  and  sinks  superimposed. 

Finally,  the  velocity  increments  given  by  the  Prandtl-Glauert 
and  Karman-Tsien  rules,  (reference  35)  are  included  in  figure  5. 

It  can  be  seen  from  figure  5  that  the  results  of  the  present  method 
differ  from  those  of  the  other  methods  by  about  the  same  amounts  as 
the  results  of  the  other  methods  differ  amongst  themselves. 

The  Prandtl-Glauert  rule  has  an  interesting  connection  with 
the  present  method.  This  rule  results  when  the  second  term  in 
equation  (34)  is  negligible  relative  to  the  first  term,  as  is  the 
case  in  the  subsonic  speed  range  for  small  perturbations.  The  tran¬ 
sonic  term  is  also  negligible,  however,  for  higher  speeds  or  larger 
perturbations  if  r  ~  2 .  As  previously  noted,  • r  — *  2  at  those  points 
on  the  airfoil  where  Ai  ~  0  but  YCa  /  0.  Near  these  (locally 

subsonic)  locations  the  Prandtl-Glauert  rule  may  therefore  be  expected 
to  apply.  On  the  other  hand,  the  maximum  velocity  increments  on  air¬ 
foils  of  the  type  illustrated  in  figure  5  turn  out  to  vary  with  Mach 
number  in  the  subsonic  range  in  a  manner  approximating  the  Karman- 
Tsien  rule.  (See  also  reference  29).  Inasmuch  as  the  KArmAn-Tsien 
rule  itself  approaches  the  Prandtl-Glauert  rule  for  small  velocity 
increments ,  the  present  method  therefore  gives  results  in  the  sub¬ 
sonic  speed  range  not  greatly  different  from  those  obtained  by  the 
Karman-Tsien  rule. 

Potential-limit  phenomenon.  -  In  computing  a  velocity  distri¬ 
bution  on  an  airfoil  at  a  fixed  free-stream  Mach  number  Mq  by 
equation  (34)  as  plotted  in  figure  3,  a  portion  of  the  airfoil  con¬ 
tour  is  put  into  correspondence  with  a  portion  of  the  lower  branch 
of  figure  3.  For  example,  on  the  Kaplan  section  at  Mq  =  0,  the 
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region  of  airfoil  0.178  <  x  <  0.823  corresponds  to  the  portion  AD 
in  figure  3,  such  that  traversing  the  stretch  0.178  <  x  <0.5  on 
the  airfoil  corresponds  to  traversing  the  stretch  AD  and  traversing 
the  stretch  0.5  <  x  <  0.823  on  the  airfoil  corresponds  to  trav¬ 
ersing  the  stretch  DA. 

As  the  Mach  number  increases,  the  same  region  of  the  airfoil 
corresponds  to  a  greater  portion  of  the  lover  branch  of  figure  3, 
the  point  A  corresponding  to  x  =  0.178  and  0.823  remaining  at  the 
origin  but  the  point  D  corresponding  to  x  =  0.5  moving  away  from 
the  origin.  For  example,  the  fixed  airfoil  location  x  =  0.5  corre¬ 
sponds  at  Mq  =  0.5  to  point  D^,  at  Mq  =  0.75  to  point  Dg,  and 

at  Mq  =  M0^j  =  0.807  to  the  turning  point  B  between  lover  and 

upper  branches  in  figure  3.  Thus,  at  a  certain  lowest  free- stream 
Mach  number  MQ  ^  a  point  x^  on  an  airfoil  comes  into  corre¬ 
spondence  with  the  turning  point  B  in  figure  3.  For  higher  but 
still  subsonic  free- stream  Mach  numbers,  the  stretch  AB  remains 
fixed  at  its  maximum  extent  but  now  corresponds  to  a  diminishing 
region  of  airfoil,  such  that  the  gap  on  the  airfoil  for  which  no 
solutions  for  the  velocity  are  yielded  by  equation  (34)  or  figure  3 
spreads  out  from  x^  as  Mc  increases  in  the  range  Mq  ^  <  Mq  <1. 

For  example,  on  the  Kaplan  section  the  forbidden  region  at  MQ  -  0.85 
is  0.35  <  x  <0.65.  Thus,  M0  ^  is  the  maximum  free-stream  Mach 

number  for  which  continuous  potential-flow  solutions  are  obtainable 
by  the  present  method.  Herein  Mc^  is  termed  nthe  potential-limit 

Mach  number,"  and  the  corresponding  flow  pattern,  "the  potential- 
limit  solution." 


Quantitatively,  the  turning  point  B  corresponds  to  the  condition 

dx/dX  =0  in  equation  (38)  or,  equivalently,  the  condition 

(du/dA)r  YC  =  const  =  0  in  eQuation  (34).  Hence,  at  the  potential 
,  a  A 

limit  (denoted  by  the  subscript  2), 


* 


2 


(46a) 


or,  by  equation  (37) 


3-r 

2-r 


(46b) 


2 


(47a) 
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21  (2-r)5 

4  /~  \  2 

(3-r) 


(47b) 


The  procedure  for  obtaining  the  potential-limit  Mach  number  for  a 
symmetric  airfoil  at  zero  lift  is  then  to  plot  the  right  side  of 
equation  (47b)  as  a  function  of  chordwise  location  x.  The  potential- 
limit  chordwise  location  x^  is  that  where  the  right  side  attains  a 

maximum  value.  This  maximum  value  yields ,  by  equation  (47b ),  a  value 
3  2 

of  |T  IK  that,  by  equations  (lib)  and  (lie),  determines  the 
potential-limit  Mach  number  M^j*  The  corresponding  velocity 
increment  is  obtained  from  equation  (47a).  The  velocity  dis¬ 

tribution  over  the  rest  of  the  airfoil  at  Mach  number  M^  is 
determined  as  usual  from  equation  (34)  or  (38)  using  figure  3. 


The  potential-limit  velocity  distributions  for  the  Kaplan  sec¬ 
tion  by  the  small-  and  finite-perturbation  integral  methods  are 
shown  in  figure  4.  The  potential- limit  Mach  numbers  by  these  two 
methods  are  0.807  and  0.843,  respectively.  Kaplan  (reference  6) 
estimates  from  the  first  three  terms  of  his  expansion  for  the  poten¬ 
tial  that  the  entire  series  would  diverge  beyond  Mq  =  0.83.  Gothert 

and  Kawalki  (reference  17)  also  locate  a  limiting  Mach  number  above 
which  their  method  does  not  converge  to  a  self-consistent  solution 
regardless  of  the  number  of  iterations  performed.  Although  the  var¬ 
ious  determinations  of  limiting  Mach  number  MQ^  and  the  corre¬ 
sponding  velocity  increment  A^,  indicated  by  the  circled  points 

in  figure  5,  do  not  agree  extremely  closely,  nevertheless,  consider¬ 
ing  the  diversity  of  the  methods  used,  the  inference  appears  war¬ 
ranted  that  a  unique,  upper-limiting,  subsonic  free- stream  Mach  num¬ 
ber  exists  beyond  which  a  continuous  potential  flow  past  a  given 
airfoil  cannot  be  derived. 


Various  properties  of  the  potential-limit  flow  pattern  can  be 
derived  by  the  present  method,  which  not  only  make  it  easy  to  under¬ 
stand  why  this  limiting  flow  pattern  should  exist  but  also  connect 
it  with  the  so-called  limiting-line  phenomenon  encountered  in  the 
hodograph  method.  (See,  for  example,  reference  36.)  In  the  first 
place,  at  the  potential  limit  the  local  Mach  number  M^  at  the 
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chordwise  location  x^  on  the  airfoil  must  he  greater  than  1,  and 
hence  a  local  supersonic  region  must  exist  in  the  flow  field.  By 
equations  (27)  and  (47a ), 


M-, 


4  ^M  / 5-r\ 
3  rc  \2-r/ 


1 


(48) 


In  the  permissible  range  (35)  for  r,  the  right  side  of  equa¬ 
tion  (48)  is  positive  down  to  at  least  M0^  =  0.63.  For  the  thin 

airfoils  contemplated  here,  MQ  ^  >  0.63  and  hence  Mj  >  1,  as 
indicated  by  the  constant  local  Mach  number  contours  in  figure  5. 


Secondly,  the  velocity  gradient  is  finite  and  suffers  a 

finite  discontinuity  at  x  =  x^  when  MQ  =  MQ  ^ •  Differentiation 
(denoted  by  a  subscript)  of  equation  (34)  with  respect  to  x  yields 


2 

(2-r)2 


r  A3 


(49) 


The  quantity  in  parentheses  on  the  left  side  of  equation  (49)  is 
zero  for  MQ  =  M Q)i>  by  equation  (47a).  The  right  side  of  equa¬ 
tion  (49)  is  zero  at  x  =  x^  for  MQ  =  M0^,  as  can  be  seen  from 
equation  (47a)  and  from  the  potential-limit  condition  that  n3/r  ^ 
in  equation  (47b)  be  a  maximum  at  x  =  x^ ;  that  is,  that  its 

derivative  by  x  is  zero  at  x^.  Thus,  x  in  equation  (49)  is 

given  by  an  indeterminate  form  at  MQ  =  M 0^,  x  =  Xj .  Application 

of  1’ Hospital’s  rule  to  this  indeterminate  form  and  use  of  equa¬ 
tion  (47)  yields  a  quadratic  equation  for  Ax.  For  the  usual  air¬ 
foil  shapes,  the  two  roots  of  this  quadratic  are  real  and  of  oppo¬ 
site  sign.  (For  two-way  symmetric  airfoils  such  as  in  fig.  5, 

the  two  roots  are  numerically  equal.)  Hence  A  undergoes,  in 

general,  a  finite  discontinuity  at  x  =  for  MQ  =  of  an 

amount  equal  to  the  difference  between  these  two  roots.  This  dis¬ 
continuity  in  velocity  gradient  is  indicated  for  the  Kaplan  section 
in  figure  4  at  x  =  0.5  and  M Q  j  =  0.807,  MQ^  =  0.843  by  the 

small-  and  finite-perturbation  methods,  respectively.  If 

(^C&)xx  =  rx  =  rxx  =  0  at  x  =  xl>  Mo  =  Mo,Z>  then  also  Ax  =  0 
and  no  discontinuity  occurs  in  Ax  at  the  potential-limit  point 
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in  this  case*  The  preceding  statements  for  Ax  also  hold  for  the 
fluid  acceleration  on  the  airfoil  inasmuch  as  the  acceleration  is 
given  by  (l+AjA^ 


Third,  by  eliminating  y  between  equations  (39)  and  (43)  and 
differentiating  the  result  by  x,  an  expression  for  the  veleoity 
gradient  (vx)s  along  a  streamline  in  the  flow  field  is  obtained 

that  is  a  linear  function  of  Ax,  (YCa)  ,  rx,  and  Ca  x,  with 

coefficients  dependent  on  A,  YCa,  r,  and  C&.  Furthermore,  the 

slope  of  a  constant  velocity,  or  constant  local  Mach  number,  line 
(yx)^,  multiplied  by  the  local  streamline  curvature  C,  equals 


(Vx)s>  Hence  Cux)s>  the  local  fluid  acceleration  along  a  stream¬ 


line,  and  (yx) 


are  associated  with  A  at  a  given  x  in  being 


finite  or  infinite  together. 


Fourth,  the  present  method  yields  a  second  flow  pattern  at  and 
above  the  potential  limit,  which  overlaps  the  first.  This  pattern 
is  given  by  the  upper  branch  of  the  cubic  curve,  figure  3;  for,  at 
the  potential  limit,  the  turning  point  B  is  just  reached  from  along 
the  lower  branch  AB  upon  traversing  the  airfoil  to  its  potential- 
limit  chordwise  location  x^ .  In  retracing  the  airfoil  in  the  oppo¬ 
site  direction,  the  upper  branch  BC  could  be  used.  For  Mq  >  MQ  j  , 
instead  of  a  single  chordwise  location  x^,  two  chordwise  locations 
xl and  xi f z  correspond  to  the  turning  point  B,  as  previously 
mentioned.  Between  the  two  lines  x  =  xi  ,1  ana  x  =  x2,2  as  "thus 

a  region  of  space  for  which  no  solution  can  be  found  by  the  present 
method.  Outside  this  region  two  overlapping  flow  patterns  are 
derivable,  corresponding  to  the  two  branches  in  figure  3.  At  the 
lines  x  =  Xj^,  x  =  XZ,2>  the  velocity  gradient  Aj^  is,  in  gen¬ 
eral,  infinite.  (The  left  factor  in  equation  (49)  is  zero  by  equa¬ 
tion  (47a),  but  the  right  side  is  not  zero,  because  of 

equation  (47b)  is  here  not  a  maximum.)  These  results  and  those  of 
the  preceding  paragraph  show  that  the  lines  x  =  ,i  and  x  =  x^, 

which  are  potential  lines  in  the  present  approximation,  are  envelopes 
of  constant-velocity,  or  constant  Mach  number,  lines.  Furthermore, 
it  is  plausible  that  the  lines  x  =  Xj  ^  and-  x  ~  xi  2  are  a^EO 

envelopes  of  Mach  lines  or  characteristics;  a  point  of  infinite 
velocity  gradient  in  a  supersonic  region  requires  a  spacing  of  Mach 
lines  closer  by  an  order  of  magnitude  than  does  a  point  of  finite 
velocity  gradient,  and  an  envelope  of  Mach  lines  satisfies  this 
condition. 
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The  preceding  properties  indicate  that  x  =  x^  ^  and 

x  =  }2  “ay  be  limiting  lines  in  the  sense  of  the  hodograph  method 

and  that  the  potential-limit  solution  is  the  boundary  solution  between 
continuous  potential-flow  patterns,  for  Mq  <  Mg  ^ ,  and  flow  patterns 

with  limiting  lines  bounding  forbidden  regions,  for  MQ  >  MQ^ .  This 

inference  is  supported  by  a  closer  examination  of  the  limiting-line 
phenomenon  as  follows: 


A  limiting  line  in  a  flow  field  is  defined  as  the  locus  of 
points  at  which  the  Jacobian  J  =  d($, V)/b(vj9)  from  the  hodograph 
(v,0) -plane  to  the  physical  ($ ,W) -plane  vanishes.  (Constant  <5 
and  ¥  define  the  potential  and  streamlines,  respectively,  in  the 
physical  plane.)  Consider  the  reciprocal  of  this  Jacobian 


„-l  d(v,0)  _  bv  bO  bv  b6 

J  “  3  (&,!!)  ~  3$ S?  SijfcW 


(50) 


The  differential  equations  of  motion  in  intrinsic  form  are  (refer¬ 
ence  35) 


b0  (M2-!)  bv 
cW "  pv  3$ 


(51a) 


S©  _  £  ^V 
3$  v 


(51b) 


Substituting  equations  (51)  into  equation  (50)  and  noting  that 
d$  =  vds  and  d¥  s  pvdn  yields 

j-1  s  pv3  J-1  =  (/-l)  (| ij  -  V2  (§§ J  (52) 

where  5v/5s  is  the  velocity  gradient  along  a  streamline  and  -50/d s 
is  the  streamline  curvature. 

The  variation  of  j"1  along  the  t  =  0,1  biconvex  airfoil, 
determined  by  equation  (52)  and  the  velocity  distributions  previ¬ 
ously  calculated,  is  shown  in  figure  6  for  various  values  of  MQ* 

Figure  6  shows  that  the  present  method  yields  finite  values  of 

in  the  range  MQ  <  The  only  discontinuity  in  behavior 

occurs  at  x  =  x^  =  0.5  at  which  point  undergoes  a  finite 
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discontinuity  for  MQ  =  M ("because,  as  previously  discussed, 
bvfbs  undergoes  here  a  finite  discontinuity  as  a  function  of  MQ) . 

It  may  he  inferred,  in  general,  that  j-*  <  0  everywhere  for 
For,  equation  (52)  can  he  solved  for  (dv/d0)s  as 

(subscript  s  denotes  values  along  a  streamline) 


-1 


(Se/ds)2 


(53) 


For  j'1  =  0,  equation  (53)  yields  the  Prandtl-Meyer  rate  of  change 
for  (bv/b9)s  in  supersonic  regions.  For  j_1>0,  equation  (53) 

yields  a  numerically  greater  rate*  Now,  in  a  local  supersonic  region  t 
embedded  in  a  subsonic  flow  field,  it  is  known  that  (Sv/d8)s  must 

be  numerically  less  than  the  Prandtl-Meyer  value  (references  37  and 
38)  because  the  sonic  boundary  reflects  Mach  waves  of  opposite  sign 
to  those  inSident  upon  it.  Flows  having  finite  fluid  accelerations, 
and  streamline  curvatures  and  containing  local  supersonic  regions 
must  therefore  have  a  finite  <  0  everywhere  (in  subsonic 
regions  j-1<  0  hy  equation  (52)).  Hence,  also  j  <0  within 
such  flow  fields.  The  preceding  discussion  may  he  regarded  as  an 
indication  of  a  simple  proof  of  a  theorem  of  Friedrichs  (reference  39) 
to  the  effect  that  if  a  continuous  potential-flow  solution  v(x,y,MQ), 
0(x,y,Mo)  varies  continuously  with  Mq  and  approaches  a  limit  as 
Mq  approaches  some  limiting  value,  then  this  limiting  solution 

itself  cannot  contain  any  points  at  which  j  =  0;  that  is,  a 
limiting  line. 


The  relatively  small  region  of  positive  values  of  j-1  indicated 
in  figure  6  for  MQ  <  MQ^  must  therefore  he  regarded  as  due  to,  and 

indicating  the  extent  of ,  the  approximations  in  the  present  method. 

A  plot  of  j~l  would  appear  to  he  a  useful  indication  for  assessing 
and  possibly  improving  the  accuracy  of  transonic  velocity  distribu¬ 
tions.  The  finite  discontinuity  in  j“-*-(M0)  at  x  =  x-^  for 

Mo  =  Mo,Z  W  wel1  exist  (in  theory).  The  present  method  thus 
illuminates  the  theorem  of  Friedrichs. 

For  M0>  M Q)i,  figure  6  shows  that  j-1  — ->»,  hence  j 0 
at  a  point  xz,i(M0)  and  at  the  symmetrically  disposed  point 
x2 ,2  =  1-XZ,1  (because  here  dv/ds — ><»,  as  previously  noted). 
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Thus,  formal  continuation  of  the  present  method  to  MQ  >  MQ^  yields 
discontinuous  solutions  containing  forbidden  regions  x^  ^  <  x  <  x^  ^ 
bounded  by  limiting  lines,  as  previously  inferred. 

The  reason  for  the  existence  of  two  values  of  X  for  a  given  X 
in  figure  3,  and  hence  two  overlapping  flow  patterns,  of  which  one 
corresponds  to  either  subsonic  or  slightly  supersonic  conditions 
and  the  other  to  highly  supersonic  conditions,  is  essentially  sim¬ 
ilar  to  the  reason  for  the  existence  of  two  possible  velocities  for 
a  given  flow  area  in  one -dimensional  flow,  as  in  a  converging- 
diverging  nozzle.  Indeed,  the  present  method,  which  makes  use  of 
the  equation  of  continuity  in  quasi -one- dimensional  or  hydraulic 
fashion,  shows  rather  directly  that  the  potential-limit  solution 
occurs  for  essentially  the  same  reason  that  a  converging- diverging 
nozzle  chokes  when  a  sufficiently  high  but  subsonic  Mach  number 
occurs  ahead  of  the  throat  (reference  29).  (The  mass-flow  density 
pv  passes  through  a  maximum  at  a  local  Mach  number  of  1.)  The 
manner  by  which  steady-flow  patterns  can  nevertheless  exist  for 
the  subsonic  free- stream  Mach  numbers  between  MQ  ^  and  unity,  in 

spite  of  the  analog  of  this  choking  phenomenon  (which  prevents  higher 
subsonic  Mach  numbers  from  being  obtained  ahead  of  the  throat  in 
nozzles)  is  considered  in  the  next  section.  In  particular,  it 
will  be  seen  that  the  purely  mathematical  extrapolation  made 
herein  leading  to  the  limiting  lines  x  =  x^  ]_  and  x  =  x^ 

has  no  physical  significance  in  connection  with  the  shock  that* 
actually  appears  at  these  Mach  numbers. 


LOWER  TRANSONIC  ASYMMETRIC  FLOW  WITH 
TERMINAL  SHOCK 

Velocity  distributions.  -  In  this  section,  asymmetric  flow 
patterns  with  a  single  terminal  shock  will  be  derived  for  free- 
stream  Mach  numbers  from  the  potential-limit  Mach  number  to  unity. 

The  analysis,  particularly  that  connected  with  the  terminal  shock, 
is  more  uncertain  in  this  range  than  in  the  neighboring  speed  ranges 
because  the  terminal  shock  is  no  longer  at  most  a  boundary  of  the 
region  of  interest, but  exists  within  it. 

The  procedure  for  obtaining  a  symmetric  solution  at  the 
potential- limit  Mach  number,  derived  in  the  preceding  section,  indi¬ 
cates  how  asymmetric  solutions  may  be  obtained  at  higher  Mach  numbers. 
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Thus,  in  determining  the  airfoil  velocity  increments  at  the 
potential- limit  free- stream  Mach  number  as  a  function  of  chordwise 
location,  the  lower  branch  of  the  curve  in  figure  3  was  traversed 
from  A  to  B,  the  potential-limit  point  B  being  reached,  say,  at  the 
midchord  location  for  the  symmetric  biconvex  airfoil.  In  going  from 
midchord  toward  the  trailing  edge  on  the  airfoil,  the  curve  was 
traced  back  from  B  to  A,  giving  the  symmetric  solution.  Because  the 
turning  point  B  at  midchord  has  just  been  attained,  however,  it  is 
possible  to  continue  around  and  have  part  of  the  upper  branch  BC 
correspond  to  the  rear  part  of  the  airfoil.  The  chordwise  location 
x  =  Xj  in  the  present  speed  range  therefore  has  the  property  of  a 

branch  point,  as  previously  explained.  This  property  maintains  the 
analogy  with  conditions  at  the  throat  of  a  converging- diverging 
nozzle. 

The  preceding  procedure  yields  an  asymmetric  solution  at  the 
potential-limit  Mach  number.  In  order  to  obtain  a  solution,  and  in 
particular  an  asymmetric  solution,  without  forbidden  regions  at 
higher  subsonic  Mach  numbers,  the  basic  cubic  equation  (34)  must  be 
modified.  The  modification  made  herein  is  to  introduce  a  nonvan¬ 
ishing,  lateral  streamline  displacement  5  at  infinity  (fig.  2). 

Such  a  displacement,  of  the  order  of  the  airfoil  thickness,  is 
clearly  necessary  for  steady-flow  patterns  to  exist  near  sonic  speed 
if  the  free- stream  flow  is  to  "get  by"  the  airfoil  in  accordance 
with  steady-state  continuity  relations.  (That  a  steady-flow  solution 
is  yielded  by  the  equations  of  motion  at  precisely  sonic  speed  has 
been  demonstrated  by  Guderley,  references  38,  40  and  41.)  The  stream¬ 
line  at  infinity  can  still  be  considered  straight  and  parallel  to  the 
x-axis,  relative  to  the  airfoil  slope,  if  its  lateral  displacement 
8  is  assumed  to  occur  gradually  over  an  infinite  distance  along  the 
streamline.  Thus,  8  will  not  be  a  function  of  airfoil  chordwise 
location  but  will  be  constant  for  the  entire  airfoil.  It  will,  how¬ 
ever,  be  a  function  of  free-stream  Mach  number.  The  introduction  of 
a  nonvanishing,  lateral  streamline  displacement  8  in  the  present 
speed  range  M.  <  M,_  7  <  1  does  not  conflict  with  the  proof  pre- 

viously  indicated  that  6  vanishes  for  subsonic-flow  conditions  and 
no  wakes  at  infinity;  the  simultaneous  introduction. of  a  shock  ter¬ 
minating  an  asymmetric  local  supersonic  region  implies  a  shock  wake 
of  slightly  decreased  total  head  extending  to  infinity.  Such  a  wake 
provides  a  mechanism  for  producing  a  nonvanishing  displacement  6 
similar  to  that  provided  by  an  ordinary  subsonic  wake  produced  by 
boundary- layer  action. 
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By  subtracting  a  suitable  displacement  function  8(MQ)  from 

Y  in  the  basic  cubic  equation  (34),  both  branches  in  figure  3  can 
now  be  correlated  continuously  with  the  airfoil  contour  in  the 
speed  range  M0^  <  Mq  <  1  to  yield  an  asymmetric  solution.  In 

order  actually  to  calculate  this  solution,  the  variation  of  the 
curvature  parameter  r  with  x  and  K,  or  its  equivalent,  must 
first  be  determined.  To  this  end,  the  assumption  is  first  made 
that  the  potential-limit  condition  of  equation  (47),  with  Y 
replaced  by  Y  -8,  will  apply  throughout  the  speed  range 
MQ  j  <  Mq  <  1  to  the  same  chordwise  location  x  =  first  found 

to  require  the  potential-limit  condition.  In  more  geometric  terms, 
if  the  basic  cubic  equation  (34)  is  viewed  as  a  surface  in  the  A, 
x,  and  MQ  space,  then  the  potential-limit  condition,  which  has 

been  regarded  as  determining  an  absolute  maximum  of  M0  on  the 
projections  of  this  surface  in  the  A,  MQ  plane,  can  therefore 

be  regarded  as  determining  the  saddle  point  of  this  surface.  Hence, 
in  the  A,  x  plane,  the  potential-limit  point  is  a  branch  point. 
Secondly,  it  is  noted  that  the  velocities  obtained  along  the  rear 
portion  of  the  airfoil,  corresponding  to  the  upper  branch  in  fig¬ 
ure  3,  and  using  the  subsonic  values  of  r  previously  determined, 
become  -much  greater  than  those  corresponding  to  a  Prandtl-Meyer 
distribution.  As  pointed  out  by  Tsien  and  Fejer  in  reference  37, 
however,  the  velocity  in  a  local  supersonic  region  should  be  less 
than  that  in  a  Prandtl-Meyer  expansion  through  the  same  angle  of 
turn,  the  difference  being  due  to  compression  waves  reflected  from 
the  sonic  boundary  and  striking  the  airfoil.  Inasmuch  as  these 
compression  waves  will  probably  be  of  small  influence  compared  to 
the  expansion  waves  emanating  from  the  airfoil  for  the  size  and  the 
shape  of  local  supersonic  regions  to  be  expected  in  the  present 
speed  range,  it  will  be  assumed  that  the  supersonic  portion  of  the 
velocity  distributions  to  be  derived  will  be  Prandtl-Meyer  dis¬ 
tributions.  (it  would  not  be  difficult  to  obtain  more  accurate 
velocity  distributions  in  local  supersonic  regions  on  the  basis  of 
the  shape  and  extent  of  the  supersonic  region  determined  by  the 
first  approximation  considered  herein.)  The  Prandtl-Meyer  velocity- 
"  distribution  formula  can  be  written  in  the  transonic  small- 
perturbation  limiting  case  as  (reference  30) 


M2-l  =  H  +  2fA  = 


3r(e-ee) 


2/3 


(54) 


where 

M-  S  -|!  =  M02-l 


(55) 
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and  0-0e  is  the  (positive)  angular  change,  in  radians,  of  slope  of 
airfoil  contour  from  the  value  0e  at  which  the  local  Mach  number 

is  1  to  the  value  9  at  which  it  is  M.  Throughout  this  section 
and  the  next,  for  simplicity  in  the  calculations,  all  T's  that 
arise,  such  as  the  of  equation  (28),  which  should  appear  in 

equation  (54),  are  replaced  by  T  =  1.2  unless  otherwise  mentioned. 
Hence,  the  calculations  are  limiting  calculations  for  Mq  =  1;  that 
is,  the  transonic  approximation  is  made  more  strictly  than  in  the 
preceding  section,  where  the  functions  Fc(Mq)  and  r^(M0)  were 
used. 


By  assuming  a  Prandtl-Meyer  distribution,  the  problem  of  cal¬ 
culating  the  supersonic  portion  of  velocity  distributions  in  a  given 
speed  range  is  reduced  to  determining  one  point  of  this  distribution 
(that  is,  velocity  and  airfoil  location)  at  each  free-stream  Mach 
number  Mq  in  the  range,  which  is  necessary  to  evaluate  the  constant 

0e  in  equation  (54).  This  one  point  will  be  chosen  as  the  velocity 
at  the  potential-limit,  or  branch-point,  location  x  =  x^  in  the 

present  speed  range.  It  is  convenient  to  carry  out  the  determination 
with  the  basic  cubic  equation  (34),  as  modified  by  inclusion  of  8, 
and  the  potential-limit  condition  of  equation  (47)  in  explicity 
similarity  form,  thus: 


where 


i  Kir  LJ 

2-r  “  3-r 

(56) 

IhJ]  (3-r) 

\2-r  J 

(57a) 

\3 

• 

(Kj  \  2 

'[&) 

(57b) 

4/(tD2/3 

1 

(58a) 

rA/(TD2/3 

(58b) 

Explicit  calculations  will  be  made  for  the  thin  symmetric  biconvex 
airfoil.  This  airfoil  contour  is  given,  for  small  thickness  ratio 
T,  by 
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!  -  2x-l  =  -  fT 


(59a) 


y  (i-e2) 


(59b) 


The  left  side  of  equation  (56)  becomes ,  by  equations  (59) 


The  potential- limit  location  for  the  symmetric  biconvex  airfoil 
is  xz  =0.5  or  |z  =  0.  Inasmuch  as  a  unique  steady-flow  pattern 

at  Mq  =  1  is  assumed,  which  is  continuous  with  the  flow  patterns 
existing  on  either  side  of  Mq  =  1,  the  velocity  at  £  =  £z  =  0 

in  particular  must  have  this  property.  In  the  next  section,  the 
1/3 

value  Lz  =  2  1  is  derived  for  the  velocity  parameter  at 
|  =  |=0  at  K  -  0  or  M0  =  1,  by  an  analysis  for  MQ  >  1. 

In  order  for  equation  (57a)  to  yield  this  value  also,  it  is  nec¬ 
essary  that  r— >2  as  K — >0  and  in  such  a  way  that 


Kl 

lim 

v  2-r 

K— >0 


3 

W = 


1.890 


(61) 


At  sonic  speed,  a  nonvanishing  streamline  displacement  at  infinity 
6  is  thereby  obtained,  given  by  equations  (57b),  (61),  and  (60) 
with  £.=  0  as 

26/t  =  0.5  (62) 

This  value  of  6,  half  the  maximum  airfoil  ordinate,  is  of  the 
right  order  of  magnitude  inasmuch  as  at  precisely  sonic  speed  a 
streamline  displacement  at  infinity  somewhat  greater  than  the  max¬ 
imum  airfoil  ordinate  is  clearly  required  by  steady- state  continuity 
considerations . 

The  velocity  at  the  potential-limit  chordwise  location  |z  =  0 
can  now  be  obtained  as  a  function  of  MQ  by  an  interpolation  pro¬ 
cedure.  The  value  of  k/( 2-r)  at  MQ  =  =  0.831  and 

!=£z=0  is  1.60  from  the  results  of  the  preceding  section.  This 


value  does  not  differ  greatly  from  that  at  MQ  =  1  given  by  equa¬ 
tion  (61).  Accordingly,  a  linear  variation  of  K/(2-r)  at  |  =  0 
as  a  function  of  K  can  be  assumed  in  this  range  that  yields  r(K) 
and  thence  L^(K)  and  &(K)  at  |  =  0,  by  equations  (57)  and 

(60).  The  values  L^(K)  at  |  =  0  enable  Prandtl-Meyer  velocity 

distributions  to  be  determined  over  the  airfoil  at  each  Mach  number 
Mq  corresponding  to  K  by  equation  (54). 

As  a  possible  refinement,  the  velocity  distribution  over  the 
front  portion  of  the  airfoil  can  be  interpolated  between  the  Prandtl 
Meyer  distribution  at  MQ  =  1  and  the  potential- limit  distribution 

at  Mq  =  Mq^  found  in  the  preceding  section.  The  interpolation  is 

made  at  each  chordwise  location  x  between  the  values  of  K/(2-r) 
given  by  equation  (56)  using  equation  (62)  and'the  values  of  L 
given  by  the  Prandtl-Meyer  distribution  at  Mq  =  1,  and  the  values 

of  K/(2-r)  at  MQ^  =  0.831  determined  in  the  preceding  section. 

The  resulting  values  of  r(x,K)  and  &(K)  determined  from  equa¬ 
tion  (57b)  then  yield  intermediate  velocity  distributions  by  equa- 
equation  (56). 

The  asymmetric-type  velocity  distributions  determined  as  just 
outlined  are  so  far  incomplete  in  two  basic  respects.  First,  the 
subsonic  portion  from  the  leading  edge  to  the  point  where  free- 
stream  velocity  is  first  reached  locally  is  lacking.  This  lack, 
as  explained  in  the  preceding  section,  is  mainly  due  to  the  inade¬ 
quacy  of  the  simple  curvature  function  (equation  (33))  in  this 
region.  For  the  purpose  of  calculating  pressure  drag  in  the  tran¬ 
sonic  range,  this  portion  is  not  needed.  (See  section  PRESSURE  EERAG 
USE  OF  MOMENTUM  INTEGRAL.)  The  second  and  more  important  omission, 
the  terminal  shock,  is  discussed  next. 

Terminal  shock.  -  The  terminal  shock  is  observed  to  form  the 
downstream  boundary  of  the  local  supersonic  region  and  to  move 
toward  the  trailing  edge  as  the  free- stream  Mach  number  increases 
toward  unity.  Assume  for  the  present  that  this  shock  is  normal 
rather  than  oblique.  This  assumption  appears  valid,  at  least  under 
the  experimental  conditions  reported  in  reference  2,  in  which  the 
measured  pressure  rises  across  the  shock  and  outside  the  region  of 
shock-boundary- layer  interaction  corresponded  fairly  closely  to 
those  of  ideal  (that  is,  Rankine-Hugoniot)  normal  shocks. 

The  velocities  across  an  ideal  normal  shock  are  related  by 
Prandtl* s  equation,  derivable  from  the  Rankine-Hugoniot  relations. 
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vxv2  =  a*2  (63) 

where  a*  is  the  speed  of  sound  at  sonic  speed  M  =  1  in  an 
isentropic  flow.  In  the  small-perturbation  transonic  case,  equa¬ 
tion  (63)  becomes,  for  the  shock  element  near  the  airfoil 
(appendix  C), 


TA^  +  TAg  =  n  (64) 

In  terms  of  the  local  Mach  numbers  M]_  3  aJ  1  +  and  Mg  s  aJi  +  jig 

upstream  and  downstream  of  the  shock,  respectively,  equation  (64) 
becomes  by  equations  (27)  and  (28b) 

dl  +  h2  =  6  (65) 


Equations  (64)  and  (65)  are  unsatisfactory  for  determining  terminal- 
shock  locations;  because  first  if  no  boundary  layer  is  assumed  and 
subsonic  values  of  A  g  are  determined  either  by  a  subsonic  rule  as 

in  reference  34,  or  equivalently  by  the  subsonic  results  of  the  pre¬ 
ceding  section,  then  the  chordwise  locations  at  which  such  subsonic 
values  of  Ag  will  "match"  the  supersonic  values  of  A-^  previously 

determined  to  satisfy  equation  (64)  or  (65)  are  within  about  20  per¬ 
cent  of  the  chord  length  of  the  trailing  edge  for  the  symmetric 
biconvex  airfoil.  This  result  would  probably  not  be  substantially 
altered  by  a  more  exact  calculation  of  viscosity-free  flow  patterns 
with  single  terminal  shock,  such  as  by  the  relaxation  method  of, 
reference  18.  Experiment  shows,  however,  (references  3  and  4)  that 
a  well -developed  shock  exists  forward  of  the  aforementioned  location. 
On  the  other  hand,  if  a  boundary  layer  is  assumed  that,  in  traversing 
the  shock,  separates  so  completely  as  to  make  the  pressure  downstream 
of  the  shock  approximately  equal  to  free- stream  pressure,  then  FAg 

(the  velocity  outside  the  boundary  layer)  drops  out  of  equation  (64). 
The  shock  locations  determined  by  TA^  =  (i,  however,  move  forward 

instead  of  rearward  as  sonic  free- stream  speed  is  approached  from 
below.  An  intermediate  situation,  in  which  Ag  could  be  determined 

from  velocity  distributions  calculated  to  allow  for  the  presence  of  a 
boundary  layer  that  has  interacted  with  the  shock,  might  yield  sat¬ 
isfactory  terminal -shock  locations  by  equation  (64).  Sue]}  a  calcula¬ 
tion  of  Ag  has  not  yet  been  attempted  because  of  the  difficulties 

of  calculating  the  shock-boundary- layer  interaction  and  its  effect 
on  the  downstream  velocity  distribution. 
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Some  understanding  of  the  rearward  movement  of  the  terminal 
shock  can  be  had,  however,  from  a  modification  of  equation  (64). 
Assume  that  the  actual  pressure  rise  across  the  normal  shock  differs 
from  the  Rankine-Hugoniot,  or  ideal  value,  by  an  amount  Ap  .  This 

difference  can  be  caused  inside  the  region  of  shock-boundary- layer 

interaction  by  friction  processes,  in  which  case  the  actual  pressure 

rise  is  less  than  the  ideal  value.  A  pressure  drop  Ap  can  also 

s 

occur  outside  the  boundary  layer,  according  to  references  2,  18, 
and  19.  This  drop  occurs  over  a  short  distance  just  downstream  of 
the  (ideal)  normal  shock  and  exists  because  of  the  necessity  for 
the  flow  to  make  a  quick  reversal  of  curvature  from  the  negative 
value  given  it  by  the  shock  to  the  positive  value  required  by  the 
airfoil  contour.  Farther  out  along  the  shock,  this  pressure  drop 
becomes  a  pressure  rise  for  the  same  reason  (reference  2).  In  any 
event,  this  type  of  pressure  change  must  physically  "join  on"  con¬ 
tinuously  with  that  in  the  boundary- layer  region.  In  view  of  the 
preceding  discussion,  a  shock  pressure-drop  coefficient  cs,  rep¬ 
resenting  the  over-all  departure  of  a  normal-shock  element  from 
Rankine-Hugoniot  behavior,  can  be  defined  in  the  small-perturbation 
transonic  case  as 


Ap. 


cs  = 


2°1V1 


2poVo 


2poVo 


=  “2V* 


2 


(66) 


If"  the  element  of  shock  front  under  consideration  is  near  the  air¬ 
foil,  as  is  the  case  here,  then  Pg  in  equation  (66)  is  the  pres¬ 
sure  coefficient  at  the  airfoil  just  downstream  of  the- shock  and 
A2  is  the  Rankine-Hugoniot  velocity  that  would  exist  just  down¬ 
stream  of  the  shock  if  cs  were  zero.  Substitution  for  from 

equation  (64)  yields 


The  additional  assumption  is  now  made  that  the  pressure  on  the 
downstream  side  of  the  shock  equals  free-stream  pressure,  or 
^2  =  0.  (in  &  conversation  with  the  author,  von  Karman  mentioned 

having  used  this  condition  some  years  ago.)  This  condition,  in 
general,  requires  separation  of  che  boundary  layer  as  it  traverses 
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the  shock.  If  cs  were  known  as  a  function  of  M*j_  (and  of  a 
suitable  Reynolds  number  at  the  shock  if  cs  represents  a  fric¬ 
tional  effect),  shock  locations  could  be  determined  from  equa¬ 
tion  (67)  with  P2  =  0.  In  the  absence  of  such  knowledge,  however, 

an  assumed  rearward  shock  movement  can  be  used  to  calculate  cs. 

Thus  assume  the  following  boundary  conditions  on  the  terminal- shock 
location:  (a)  At  the  potential-limit  Mach  number  MQ  the  ter¬ 
minal  shock  occurs  at  the  potential-limit  chordwise  location  x^ 

(x^  *s  0.5  for  the  biconvex  airfoil);  and  (b)  At  a  free-stream  Mach 
number  of  unity,  the  terminal  shock  is  at  the  trailing  edge,  x  =  1. 

These  boundary  conditions  are,  of  course,  approximate.  The 
terminal  shock  appears  to  be  formed  by  the  ever-present  random  pres¬ 
sure  disturbances  that  are  continually  propagating  through  the  flow 
field  at  locally  sonic  speed.  Some  of  these  disturbances,  traveling 
upstream  (reference  l),  encounter  the  locally  supersonic  region  and 
tend  to  pile  up;  that  is,  the  terminal  shock  appears  to  form  as  an 
envelope  of  characteristics  in  the  x,y,t  space  (t  is  time)  rather 
than  as  an  envelope  of  characteristics  in  the  x,y  plane  (limiting 
line).  It  may  be  necessary,  however,  for  the  symmetric-type  locally 
supersonic  region  to  reach  a  finite  size  before  the  unsteady  dis¬ 
turbances  can  coalesce  to  form  a  steady  terminal  shock  (reference  4). 
The  boundary  condition  (a)  idealizes  this  process  by  assuming  no 
shock  and  a  symmetric  solution  up  to  the  potential-limit  Mach  num¬ 
ber,  thence  a  sudden  change-over  to  an  asymmetric  solution  with 
steady  terminal  shock.  Presumably,  the  disturbances  that  form  the 
terminal  shock  have  a  tendency  to  pile  up  as  far  upstream  as  possible. 
The  terminal  shock  could  hardly  form  upstream  of  the  potential-limit 
location,  however,  because  of  the  "choked"  condition  at  this  location. 
(See  preceding  section.  In  this  connection  compare  the  experimental 
data  and  the  semiempirical  analysis  of  reference  24.)  As  regards 
boundary  condition  (b),  there  are  no  experimental  data  on  terminal- 
shock  location  and  structure  at  Mach  number  unity.  In  the  super¬ 
sonic  speed  range,  this  shock  becomes  the  oblique  shock,  which  does 
start  close  to  the  trailing  edge.  The  schlieren  photographs  of 
Liepmann  (reference  3)  suggest  that  the  terminal  shock  will  reach 
the  vicinity  of  the  trailing  edge  at  sonic  speed  but  will  also  be 
curved  in  the  direction  ultimately  required  in  the  completely 
supersonic  regime. 

Boundary  condition  (a)  and  equation  (67),  with  P2  =  0,  now 
yield  for  the  T  =0.1  biconvex  airfoil,  at  M0  =  MQ  l  =  0.831 
and  =  0.389  (values  shown  in  fig.  5),  the  shock  pressure- 

drop  coefficient  cs  =  0.26.  At  Mq  =  1  and  x  =  1,  A^^  -  o.47 
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(from  the  following  section  and  fig.  7)  giving  cs  =  0.94.  The 
assumption  of  a  linear  variation  of  with  1 1  between  the  pre¬ 
ceding  end  values  (the  extreme  values  of  do  not  differ  greatly) 

determines,  from  the  velocity  distributions  previously  calculated, 
the  variation  of  terminal- shock  location  xg  with  Ifech  number  Mq 

(fig.  8).  Typical  velocity  distributions  as  thus  completed  by 
location  of  the  terminal  shock  are  shown  in  figure  7.  (The  velocity 
increment  A  is  related  to  the  pressure  coefficient  P  by 
P  =  -2A.  By  the  assumption  P2  =  0,  A  is  therefore  indicated  as 

zero  behind  the  shock.  This  A  refers  to  conditions  outside  the 
boundary  layer.) 

Other  shock  conditions  have  been  advanced.  Tsien  and  Fejer 
(reference  37)  regard  the  shock  as  an  oblique  shock  producing  max¬ 
imum  outward  deflection  of  the  streamline  outside  the  boundary 
layer.  Theodorsen  (reference  42)  regards  the  shock  as  the  oblique 
shock  that  produces  locally  sonic  speed  on  its  downstream  side. 

In  the  small-perturbation  transonic  case,  these  shock  conditions 
lead  to  equations  similar  in  form  to  equation  (64):  namely 
(appendix  C), 


Maximum  deflection: 

TA1/3 

+  ta2  =  2/3 

(68) 

Croccofs  point: 

r  a.j/7 

+ ta2  =  4/7 

(69) 

Sonic  point: 

ta2  =  1/2 

P- 

(70) 

An  additional  shock  condition  has  been  included  in  equation  (69). 

This  condition  is  the  oblique  shock  for  which  the  streamline  cur¬ 
vature  on  the  downstream  side  is  zero  (if  the  oblique  shock  has 
finite  curvature,  reference  43).  It  is  defined  by  Crocco's  point 
on  the  shock  polar  ( appendix  C ) . 

From  the  similarity  between  the  preceding  equations  and  equa¬ 
tion  (64),  essentially  the  same  objections  as  regards  ability  to 
predict  shock  location  that  have  been  made  to  equation  (64)  can  be 
made  to  the  preceding  equations.  In  addition,  as  has  been  mentioned, 
measurements  indicate  the  terminal  shock  to  be  normal  under  the 
experimental  conditions  of  reference  2.  Making  the  same  mod¬ 
ifications,  however,  in  equations  (68)  and  (69)  as  were  made  in 
equation  (64),  that  is,  introducing  a  shock  pressure-drop  coefficient 
cg  and  assuming  Pg  =0,  yields,  with  the  same  shock  boundary  con¬ 
ditions  as  before,  practically  the  same  variation  of  shock  location 
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with  Mach  number  M0  but  different  ranges  of  variation  of  cs. 

These  conditions  are  listed  in  table  I.  (Equation  (70)  does  not 
yield  a  solution  for  shock  location  inasmuch  as,  unlike  the  previous 
equations,  it  does  not  relate  conditions  across  the  shock  but  speci¬ 
fies  only  the  velocity  on  its  downstream  side.  Equation  (70)  there¬ 
fore  depends  in  an  essential  way  on  the  knowledge  of  the  velocity 
distribution  downstream  of  the  shock.) 


TABLE  I 


Terminal- shock  condition 

Shock-pressure-drop 
coefficient,  cg 

M  7  =  0.831 

v,  t 

M0  =  1 

Normal 

0.26 

0.94 

Maximum  deflection 

-  .09 

.31 

Crocco's  point 

-  .18 

.14 

Sonic  point 

-  .26 

0 

The  pressure  measurements  of  reference  2  indicate  values  of 
cs  of  about  0.1  just  outside  the  boundary  layer  and  of  about  0.3 

at  the  surface  of  the  airfoil  at  an  effective  free- stream  Mach  number 
probably  much  closer  to  the  potential-limit  Mach  number  than  to 
unity.  The  surface-pressure  and  wake  survey  measurements  of  refer¬ 
ence  44  can  be  interpreted  as  indicating  values  of  cg  up  to  about 

0.5  in  the  boundary- layer  region.  These  estimates  of  cs  would 
appear  to  favor  the  normal-shock  condition  (table  I),  particularly 
inasmuch  as  the  negative  values  calculated  in  table  I  are  probably 
unrealistic.  The  actual  shock-boundary- layer  structure  near  the 
airfoil  is,  of  course,  considerably  more  complicated  than  that  of 
a  single  straight  shock,  whether  normal  or  oblique.  The  introduction 
of  a  coefficient  such  as  cs  to  represent  an  over-all  departure 
from  ideal  single- shock  behavior,  however,  may  possibly  be  useful 
practically. 

The  reason  for  the  rearward  movement  of  the  terminal  shock  with 
increase  of  Mach  number  is  not  necessarily  indicated  by  the  pre¬ 
ceding  correlation  between  terminal- shock  location  and  pressure-drop 
coefficient  cs.  It  would  seem  more  likely  that  a  given  terminal- 
shock  location  resulting  from  other,  more  powerful,  dynamical 
requirements  should  cause  a  particular  value  of  cs  to  exist, 
rather  than  conversely.  These  other  dynamical  requirements  may  be 
connected  with  the  necessity  for  the  streamline  displacement  at 
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infinity  6(MQ)  to  increase  with  Mach  number  in  the  present  speed 

range.  An  increasing  6  could  be  produced  by  a  more  extensive 
and  intensive  shock-boundary- layer  wake,  which  in  turn  would  require 
a  longer  and  stronger  terminal  shock,  such  as  would  occur  at  the 
higher  local  Mach  numbers  farther  to  the  rear  of  the  airfoil. 

It  need  hardly  be  emphasized  that  the  calculations  of  this 
report,  particularly  those  in  this  section,  are  of  a  provisional 
and  exploratory  character.  At  the  present  stage,  not  much  more  has 
been  rigorously  demonstrated  than  that  the  methods  used  herein  can 
yield  results  of  the  right  order  of  magnitude. 

Finally,  it  is  of  interest  to  note  that  the  condition  used  in 
this  section,  that  shock  back  pressure  shall  equal  free- stream  pres¬ 
sure,  can  be  regarded  as  a  special  case  of  a  more  general  condition j 
namely,  that  the  shock  back  pressure  shall  be  as  close  as  it  can  get 
to  free-stream  pressure.  This  more  general  condition  is  discussed 
in  appendix  D. 


UPPER  TRANSONIC  FLOW  WITH  DETACHED  SHOCK 


A^lysis.  -  The  upper  transonic  or  detached  shock  regime  extends 

from  free-stream  Mach  number  1  to  that  supersonic  Mach  number  M„  0 

o  j  a 

at  which  the  shock  ahead  of  the  airfoil  becomes  attached  to  the 
leading  edge,  assumed  sharp-edged.  In  the  transonic  small- 
perturbation  case,  MQ>a  is  given  by  (reference  30  and  appendix  C), 


M 


o.a 


'-1 


(eTr) 


273  7273 


=  1.890 


(71) 


where  0^  is  the  semivertex  angle  at  the  leading  edge  (fig.  9). 


Calculation  will  now  be  made  of  the  movement  of  the  head  shock 
(fig.  1)  and  of  the  supersonic  part  of  the  airfoil  velocity  dis¬ 
tribution  as  a  function  of  MQ.  The  terminal  shock  is  assumed  to 

start  from  the  trailing  edge  (see  preceding  section)  so  that  it  need 
not  be  considered  in  the  analysis.  As  in  the  preceding  section,  the 
supersonic  part  of  the  velocity  distribution  on  the  airfoil  is 
assumed  to  be  a  Prandtl-Meyer  distribution  given  by  equatipn  (54). 
This  assumption  increases  in  validity  with  the  Mach  number.  At  a 
Mach  number  of  1,  the  flow  pattern  derived  by  Guderley  (refer¬ 
ence  40)  indicates  the  approximate  validity  of  such  an  assumption, 
except  in  a  small  region  near  the  beginning  of  the  supersonic 
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velocity  distribution  that  is  relatively  more  substantially  influ¬ 
enced  by  compression  vaves  from  the  sonic  boundary.  The  preceding 
assumption  is  therefore  probably  a  much  better  one  in  the  present 
speed  range  than  in  that  of  the  preceding  section. 

As  in  the  preceding  section,  by  assuming  a  Prandtl-Meyer  dis¬ 
tribution  the  problem  of  calculating  the  supersonic  part  of  velocity 
distributions  in  the  upper  transonic  range  is  again  reduced  to  deter¬ 
mining  one  point  of  this  distribution  at  each  free-stream  Mach  num- 
ber  Mq  in  the  range,  which  is  necessary  to  evaluate  the  constant 

0e  in  equation  (54).  This  calculation  is  performed  along  the  same 

lines  as  in  the  preceding  section.  The  continuity  condition  is  first 
set  up  for  the  flow  across  a  section  AB  in  the  y-direction  (fig.  9), 
where  A  is  a  point  on  the  shock,  to  be  specified  presently,  at  a 
distance  yw  from  the  x-axis  and  B  is  the  laterally  opposite  point 

on  the  airfoil  at  the  same  chordwise  location  xw .  Thus 


r^yw 

Y  =  (pv-l)dy 

J  Y 


(72) 


The  irrotationality  condition  (equation  (5))  applies  also  in  the 
present  speed  range,  inasmuch  as  fluid  rotation  introduced  by  var¬ 
iable  entropy  increase  across  a  curved  shock  in  the  small-perturbation 
transonic  case,  being  of  the  order  of  the  third  power  of  the  Mach  num¬ 
ber  increment  across  the  shock,  is  negligible.  Substituting  equa¬ 
tion  (5)  into  equation  (72)  and  making  the  small-perturbation  sub¬ 
stitution  of  equation  (10)  for  the  density  yields 


Y 


(73) 


where  vw  is  the  velocity  increment  on  the  downstream  side  of  the 
shock  at  A  and  A  is  the  velocity  increment  at  B. 

A  streamline-curvature  function  is  now  assumed  as 


C 


=  C, 


(74) 
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where  s  is  a  parameter  to  he  determined  as  a  function  of  and 

K.  Equation  (74)  satisfies  the  boundary  condition  at  the  airfoil 

and  in  addition  states  that  the  curvature  is  zero  at  point  A.  It  is 

desirable  that  the  streamline  curvature  should  be  zero  at  point  A 

because  if  it  were  not,  then  the  determination  of  the  curvature  at 

point  A  would  greatly  complicate  the  problem.  Now  a  point  A  exists 

on  a  curved  shock  at  which  the  streamline  curvature  Just  downstream 

of  the  shock  is  zero.  It  is  called  Crocco's  point  (references  38,  40, 

41,  45,  and  46)  and  is  defined  by  the  condition 

(dv/§0)  i  ,  ,  .  =  oo  (appendix  C).  Instead  of  choosing 

'  along  streamline  ■ 

point  A  as  Crocco's  point,  however,  it  will  be  chosen  as  the  sonic 
point,  that  is,  as  that  point  at  which  the  local  Mach  number  just 
downstream  of  the  shock  is  unity.  The  reason  for  this  choice  is 
the  following:  Because  of  the  additional  degree  of  freedom  pre¬ 
sented  by  the  position  of  the  shock  relative  to  the  airfoil,  that  is, 
determination  of  the  distance  yw,  one  more  condition  in  addition 

to  those  already  given  is  necessary* to  make  the  problem  determinate. 

If  point  A  is  the  sonic  point,  a  very  simple  such  condition  is 
obtained,  because  at  the  sonic  point,  (^d/^v)  ^  =  q  (appen¬ 

dix  C).  Hence,  to  the  first  order  in  the  change  of  velocity  between  A 
and  B  (that  is,  in  the  small-perturbation  case)  the  stream  direc¬ 
tion  Qy  at  A  will  equal  the  airfoil  inclination  £?a  at  B.  The 

violation  of  equation  (74),  because  the  curvature  will  not  be  zero 
at  A,  is  probably  not  serious  inasmuch  as  it  is  known  that  the  sonic 
point  and  the  point  of  maximum  streamline  deflection  are  quite  close 
to  each  other  on  the  shock  polar  and  Crocco's  point  is  between  these 
two  points  (appendix  C).  Hence  the  streamline  curvature  should  be 
close  to  zero  at  the  sonic  point  A  in  the  physical  plane.  It  is 
noted  that  choice  of  a  point  with  a  particular  property  on  the  head 
shock  to  determine  the  path  of  integration  in  equation  (72)  makes 
the  chordwise  location  x^  a  function  of  the  Mach  number  parameter  K 

in  the  analysis.  Hence  the  parameter  s  in  equation  (74)  becomes  a 
function  only  of  K. 

Substitution  of  equation  (74)  into  equation  (73)  and  integration 
yields 


=0w=' 


(m>w+r»w‘ 


(l-s) 


<A-v* 


TTiT  (A-V 
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in  which  the  aforementioned  condition  6a  =  Gv  is  indicated  in  the 
subscript  on  the  left  side. 

The  distance  yw  of  the  shock  point  A  from  the  airfoil  point  B 

is  obtained  by  integration  of  the  irrotationality  condition  (equa¬ 
tion  (5))  with  the  curvature  function  (equation  (74))  as 

A~vv  ,  , 

y*  ■  Y  -  -  ca(i-s) 

In  order  to  satisfy  the  condition  of  zero  curvature  at  V  =  by 

equation  (74)  and  for  convergence  of  the  integrals  leading  to  equa¬ 
tions  (75)  and  (76) ,  it  is  necessary  to  restrict  the  range  of  values 
for  the  parameter  s  to 


0  <  s  <  1  (77) 

The  distance  d  of  the  normal  element  E  of  the  shock  from  the 
leading  edge  0  of  the  airfoil  can  be  expressed  as  follows:  The 
angle  of  inclination  o  of  the  tangent  AC  to  the  shock  at  the  sonic 
point  A  is  given  by  (appendix  C) 


cos  O 


(78) 


This  inclination  is  somewhat  steeper  than  that  for  a  Mach  wave,  for 

which  cos  a  =  as  is  to  be  expected.  The  distance  CD  is  thus 

given  by 

CD  =  yw  cot  a  z  yw  cos  a  (79) 


Introducing  a  factor  A  to  allow  for  the  reduction  CE  due  to  cur¬ 
vature  of  the  shock  between  A  and  E  therefore  gives,  using 
equation  (76), 


d+xw  =  A  l[n/ 2  yw  =  — 


p/2  (A-Vw) 


CA(l-s) 


(80) 


The  preceding  equations  will  now  be  applied  to  the  thin  biconvex  air¬ 
foil.  For  this  contour,  conditions  at  the  chordwise  location  xw 

analyzed  here  are  given  by  equation  (59)  as 
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*-Ul-rT 


2  yw 

(YCa)  -  2t  - 

0a  =  0v  (82) 

It  is  convenient  to  proceed  with  the  equations  expressed  explicitly 
in  transonic- similarity  form.  At  the  sonic  point  behind  a  shock 
(appendix  C), 


F0v  1 


(83a) 


(83b) 


Substituting  equation  (81)  into  (80),  equation  (82)  into  (75),  noting 
that  the  second  term  on  the  right  side  in  equation  (75)  vanishes 
because  of  equation  (83a),  and  expressing  equations  (75),  (81),  and 
(80)  in  similarity  form  with  the  aid  of  equations  (83a)  and  (83b) 
and  the  definitions 

J/2 (tT)2''3  .  -K/2  (84) 

L  -  ^  A/t2^3  (*, 

yields,  respectively, 


(YCa)0„  =  0. 


K-i3  K-,2 


V  % 

_  =  2  -  —  =  —  (L^+Ki) 


(Lv+%) 


=  r  (l  -  l  K  3/2 


d  =  iiBj  (v%)  -  i  1  k. 


The  parameter  s  in  equation  (86)  is  next • determined  as  a  function 
of  Kj_.  To  begin  with,  as  K-^ — )0,  s  must— >1,  and  in  such  a  way 
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that  Kq2/(l-s) — ^constant,  for  if  not,  then  at  =  0,  Ly  <  0 

(by  equation  (86)),  which  is  incorrect  in  view  of  the  direction  of 
curvature  of  the  flow.  It  will  be  assumed  that  the  variation  of  s 
with  K]_  is  given,  not  only  near  =  0,  but  throughout  the  entire 
upper  transonic  range  0  <  <  3/2  by 

K  2 

— =  B  (a  constant)  (89) 

1-s 


The  value  of  this  constant  can  be  determined  as  follows:  Equa¬ 
tion  (86)  is  similar  in  structure  to  the  corresponding  subsonic 
equation  (56).  It  was  assumed  in  the  preceding  section  that  equa¬ 
tion  (56)  should  comply  with  the  potential-limit  condition 


(90) 


at  x  =  x^  =  0.5  throughout  the  lower  transonic  range.  The  sub¬ 
script  x  in  equation  (90)  indicates  that  x  is  held  constant  in 
the  differentiation.  Now,  if  a  unique  continuous  solution  exists  at 
MQ  =  1,  then  equation  (86)  should  become  equivalent  to  equation  (56) 

at  MQ  =  1  and  hence  should  obey  the  potential- limit  condition 
analogous  to  equation  (90)  at  MQ  =  1  or  K]_  =  0  and  x  =  0.5.  An 

analogous  condition  to  equation  (90),  rather  than  the  same  condition, 
is  necessary  because  here  x  =  Xy  is  a  function  of  and  so  can¬ 

not  be  held  constant  in  differentiating  L.  However,  inasmuch  as 


^w  (*L  \ 

<%  "  \SkI  )  + 


dxy 

dKq 


(91) 


and  the  second  term  on  the  right  side  is  finite  and,  in  fact,  zero  at 
K-^  =  0  by  equation  (87),  it  is  seen  that  the  analogous  potential- 

limit  condition  to  be  applied  to  equation  (86)  is 


=  0 


%  =  0 


(92) 
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Tlie  chordwise  location  at  which  equation  (92)  is  applied  is,  at 
Kl  =  0,  by  equation  (87),  x^  =  0.5,  which  joins  on  continuously 

with  the  branch-point  location  for  the  biconvex  airfoil  in  the  lower 
transonic  speed  range.  On  an  airfoil  for  which  the  =  0  chord- 

CL 

Wise  location  differs  from  the  branch-point  chordwise  location  x^ 

previously  derived,  interpolation  between  these  two  locations  as  a 
function  of  K  would  be  necessary  in  order  to  obtain  the  branch¬ 
point  locations  in  the  lower  transonic  asymmetric-flow  range.  Apply¬ 
ing  the  condition  of  equation  (92)  to  equation  (86)  results  in,  by 
equation  (89), 


B=^L  .  2 
2 


Inserting  this  value  into  equation  (86)  yields 


^,1  ~ 

and,  by  equation  (93), 


>—  T 


1/3 

=  2  '  .  =  1.260 


B  =  =  2.380 


P73 


(93) 


(94) 


(95) 


The  remaining  constant  A,  occurring  in  equation  (88),  is 
evaluated  by  the  boundary  condition  that  the  shock  distance  d— >  0 
as  where  K-^a  corresponds  to  the  Mach  number  for 

attachment  of  the  shock  and  is  given  by  equations  (71),  (59),  and 
(84)  as  K-^a  =  3/2.  With  these  values  for  d  and  K^_  and  by 

equations  (89)  and  (95),  equation  (88)  yields 

(^+3/2)  A  =  0.126  (96) 

On  the  other  hand,  equations  (89),  (95),  and  (86)  yield  a  cubic  for 
Ly+3/2  at  Ki=3/2,  with  the  solution 

1^+1  =  0.1318  (97) 

Hence 


A  =  0.955 


(98) 
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The  constant  A  is  of  the  expected  order  of  magnitude,  namely,  some- 
vhat  less  than  unity,  although  perhaps  too  close  to  unity  to  indicate 
a  high  degree  of  quantitative  accuracy. 

The  other  boundary  condition  on  shock  distance,  d — ) oo  as 
Kl — >0,  is  already  satisfied  by  equation  (88),  which  gives  a 

limiting  variation,  using  equation  (89)  and  the  preceding  numerical 
values , 


lim  d— > 

Kj- *0 


ABL- 


4K- 


TJz 


0.716 


(99) 


The  variation  of  s  with  K^,  given  by  equation  (89),  and  the 

various  numerical  constants  just  derived  enable  such  quantities  as 
the  longitudinal  shock  distance  d,  the  lateral  shock  distance 
yw,  the  chordwise  location  Xy,  and  the  velocity  at  chordwise 

location  x^,  on  the  airfoil  to  be  derived  as  functions  of  K-^  from 
equations  (76)  and  (86)  to  (88).  From  LW(K-^),  x^K^),  which  con¬ 

stitute  the  determination  of  a  known  point  of  the  velocity  distribution 
at  each  K^,  the  supersonic  part  of  the  velocity  distribution  on  the 
symmetric  biconvex  airfoil  is  obtained  by  the  Prandtl-Meyer 
formula  (54). 


The  various  shock  distances  calculated  as  indicated  are  shown  in 
figure  10.  Some  velocity  distributions  are  shown  in  figure  7.  The 
numerical  values  in  figures  7  and  10  are  plotted  for  a  biconvex  air¬ 
foil  of  =  0.1.  These  values,  however,  were  actually  derived  from 
the  similarity  parameters  calculated  by  the  equations  of  this  section, 
which  yield  only  limiting  values  for  j  =  0,  MQ  =  1.  As  an  indi¬ 
cation  of  the  error  involved  in  applying  the  results  of  figures  7  and 
10  to  biconvex  airfoils  of  finite  thickness  ratio,  figure  10  indicates 
the  Mach  number  for  shock  attachment  to  a  t  =  0.1  biconvex  airfoil 
as  1.31;  whereas  the  exact  oblique- shock  relations  show  it  to  be  1.48. 
Practical  application  of  the  results  of  this  section  therefore  requires 
further  investigation  in  the  direction  of  the  counterpart  of  the 
present  calculation  for  T  and  p  finitely  different  from  zero. 
(Calculations  by  integral  methods  of  detached- shock  configurations 
relative  to  thick  bodies  have  recently  been  made  by  W.  E.  Moeckel, 
reference  47.) 

Some  comparisons  with  experiment.  -  Figure  7  has  been  cross - 
plotted  in  figure  11  to  show  the  velocity  increments  A  at  the 
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x  =  0.5  and  x  =  0.8  chordwise  locations  through  the  entire  tran¬ 
sonic  speed  range.  From  these  velocity  increments,  local  Mach  num¬ 
bers  M  have  been  calculated  by  equation  (27),  assuming  TM  =  T=  1.2 

for  over-all  consistency  with  the  determination  of  A.  These  Mach 
numbers  are  shown  as  the  solid  lines  in  figure  12.  The  relative  con¬ 
stancy  of  M  at  x  =  0.5  in  the  transonic  range  may  be  especially 
noted,  (in  the  low  subsonic  and  high  supersonic  ranges,  M  should 
approximately  equal  M0.)  This  phenomenon  is  of  particular  interest 
in  view  of  the  rather  independent  derivations  of  the  maximum  local 
Mach  number  at  the  potential-limit  Mach  number  MQ  j  and  at  sonic 

speed  Mq  =  1.  The  potential-limit  derivation  depended  on  the 

incompressible  velocity  distribution  for  the  parameter  r;  whereas 
the  sonic-speed  derivation  depended  on  a  property  at  the  sonic  point 
of  the  shock  polar  (S0/Sl>)y  =  0.  The  relatively  constant  value  of 

M  at  x  =  0.5  throughout  the  transonic  speed  range  is  plausible 
from  the  expected  tendency  of  the  flow  to  maintain  a  maximum  average 
value  of  pv  ~  1  along  the  line  x  =  x^ ,  when  both  subsonic  and 

supersonic  velocities  occur  along  this  line. 

The  discontinuities  in  the  theoretical  curves  for  x  =  0.8  at 
MQ  =  0.83  and  at  Mq  =  0.94  correspond  to  the  assumed  sudden 

appearance  of  the  terminal  shock  at  the  potential-limit  Mach  number 
(and  the  condition  Pg  =  0)  and  to  the  shock  passing  through  the 

x  =  0.8  location,  respectively,  (fig.  8). 

Pressure  distributions  in  the  lower  transonic  range  on  6-  and 
12-percent  thick  symmetric  biconvex  airfoils  of  3-inch  chord  spanning 
a  two-dimensional  tunnel  20  inches  high  and  2  inches  wide  were  meas¬ 
ured  in  references  3  and  4.  The  Eeynolds  number  based  on  airfoil 
chord  was  of  the  order  of  one  million.  The  results  were  not  cor¬ 
rected  for  tunnel-wall  effect.  Terminal- shock  locations  plotted  in 
figure  8  were  estimated  herein  from  the  given  pressure  distributions 
and  schlieren  photographs.  The  local  Mach  numbers  given  in  refer¬ 
ence  4  were  presumably  obtained  by  applying  the  finite-perturbation 
formulas  (equations  (6)  and  (26))  to  the  pressure  differences  actually 
measured.  Accordingly,  these  formulas  were  used  in  obtaining  velocity 
increments  and  local  Mach  numbers  from  the  experimental  results  given. 
All  velocity  increments  A,  local  Mach  numbers  M,  and  free-stream 
Mach  numbers  MQ  were  reduced  from  T=  0.06,  0.12  to  T  =  0.10  by 

the  transonic  similarity  rules  of  equations  (24),  (25),  and  (31). 

The  experimental  results  as  thus  reduced  are  shown  in  figures  8,  11, 
and  12. 
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Several  points  uncorrected  by  the  transonic  similarity  rules 
are  included  in  figures  11  and  12.  It  is  seen  that  these  rules  do 
substantially  correlate  the  different  thickness-ratio  results.  The 
residual  discrepancy  in  the  results  for  x  =  0.5  after  correlation 
may  be  due  to  the  finite  difference  between  Mq  and  1,  to  the 
finite-perturbation  velocities  involved,  to  differential  tunnel-wall 
effects,  or  to  differential  Reynolds  number  effects.  Also  included 
in  figures  11  and  12  are  several  points  similarly  reduced  from  the 
tests  described  in  reference  48.  These  tests  were  made  on  a  6-percent 
thick,  4-inch  chord  symmetric  biconvex  airfoil  in  a  two-dimensional 
tunnel  18  inches  high  and  4  inches  wide.  Reynolds  number  was  of  the 
order  of  one  million.  The  discrepancy  between  the  highest  speed 
points  of  reference  48  and  the  lowest  speed  points  of  references  3 
and  4  is  possibly  due  to  larger  tunnel-wall  effects  and  the  somewhat 
higher  Reynolds  number  in  the  tests  of  reference  48. 

Agreement  of  the  present  theory  with  the  experimental  results 
is  considered  to  be  good  at  least  qualitatively.  The  theoretical 
calculations  are,  as  previously  noted,  limiting  calculations  for 
T  =  0  and  M0  =  1.  As  an  indication  of  the  necessity  of  allowing 
for  finite  departures  from  at  least  the  transonic  limit  Mq  =  1  in 
order  to  attain  closer  agreement  with  experiment,  several  alternative 
theoretical  local  Mach  number  curves  are  included  in  figure  12.  In 
one,  the  function  r^(Mo)  of  equation  (28a)  was  used  in  equation  (27) 

to  obtain  local  Mach  numbers  from  the  theoretical  velocity  increments 
of  figure  11.  In  the  other,  the  exact  equation  (26)  was  used.  It 
is  seen  that  the  differences  between  these  alternative  curves  is  of 
the  same  order  as  the  difference  between  theory  and  experiment.  The 
experimental  local  Mach  numbers  appear,  nevertheless,  to  be  definitely 
lower  and  higher  than  the  theoretical  values  at  x  =  0.5  and  x  =  0.8, 
respectively.  It  remains  to  be  seen  whether  experiments  at  higher 
Reynolds  number  would  not  better  the  agreement  with  the  (essentially 
infinite  Reynolds  number)  theory.  The  changes  with  increase  of 
Reynolds  number  may  be  expected  to  be  in  the  direction  for  better 
agreement. 


PRESSURE  DRAG;  USE  OF  MOMENTUM  INTEGRAL 

The  flow  patterns  obtained  in  the  preceding  sections,  together 
with  the  momentum  equation  in  integral  form,  yield  an  estimate  of 
the  pressure  drag  of  the  symmetric  biconvex  airfoil  in  the  transonic 
speed  range.  The  pressure  drag  of  an  airfoil  is  defined  as  that 
portion  of  the  total  drag  that  can  be  obtained  by  an  integration 
around  the  airfoil  of  the  component  in  the  free-stream  (x-)  direc¬ 
tion  of  the  normal  force  or  pressure  on  the  airfoil.  The  remaining 
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portion  of  the  total  airfoil  drag,  the  skin  friction,  due  to  the 
tangential  forces  on  the  airfoil  surface  produced  by  viscous  action 
in  the  boundary  layer,  will  not  be  considered  here.  The  skin  fric¬ 
tion,  for  example,  on  a  10-percent  thick  airfoil,  is  known  to  be 
small  relative  to  the  pressure  drag  in  the  transonic  range. 

In  the  range  of  continuous  potential  flow,  that  is,  from  Mach 
number  zero  to  the  potential- limit  Mach  number  MQ  j ,  the  pres¬ 
sure  drag  is  zero,  corresponding  to  the  fact  that  the  completely 
subsonic  potential-flow  field  bounding  the  airfoil  provides  no 
mechanism  for  transporting  x-direction  momentum  away  from  the  air¬ 
foil.  In  the  range  of  transonic  flow  with  shock,  however,  the  shock 
structure  provides  such  a  mechanism,  both  in  the  direct  shock  wake 
of  decreased  total  head  and  in  the  additional  wake  due  to  boundary- 
layer  thickening  in  traversing  the  terminal  shock.  In  this  range, 
a  nonvanishing  pressure  drag  is  therefore  obtained  that  could  be 
directly  calculated  by  integration  over  the  airfoil  of  the  x- component 
of  the  pressure.  This  direct  integration  can  be  carried  out  here  over 
the  rearward  part  of  the  symmetric  biconvex  airfoil  for  which  velocity 
distributions  have  been  previously  derived.  The  contribution  of  the 
remaining  forward  part  of  the  airfoil  to  the  pressure  drag  can  be 
obtained  from  knowledge  of  the  conditions  in  the  flow  field  by  use  of 
the  momentum  integral . 


Integration  of  the  x- component  of  the  pressure  between  two  chord- 
wise  locations  xa  and  x^  yields  for  the  pressure-drag  coefficient 

cd,r  this  Part  of  airfoil  (upper  and  lower  surface) 


cd,r  =  J 


2 

2povo  c 


=  2 


P  tan  0a  dx 


=  -4 


*b 


A0adx 


(100) 


in  the  small-perturbation  case.  For  the  symmetric  biconvex  airfoil 
equation  (100)  becomes,  by  equation  (59a), 


(101) 


Substituting  the  transonic  small-perturbation  limiting  form  (54)  of 
the  Prandtl-Meyer  velocity  distribution  into  equation  (101)  yields 
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The  resulting  pressure  drag  on  the  symmetric  biconvex  airfoil 
is 

cd  ==  cd,f  +  cd,r  (105) 

and  is  shown  in  figure  13  as  a  function  of  free- stream  Mach  number. 
The  finite  value  of  c^  =  0.015  at  the  potential-limit  Mach  num¬ 
ber  MQ  =  0.83  results  from  placing  the  terminal  shock  at  the 

potential-limit  chordwise  location  at  this  Mach  number  (see  pre¬ 
ceding  section) .  As  the  terminal  shock  moves  toward  the  trailing 
edge,  the  drag  rapidly  rises  to  a  maximum  value  of  c^  =  0.098 

at  sonic  speed  M^  =  1.  There  are  several  points  of  interest  about 

this  value.  First,  it  is  of  the  same  order  of  magnitude  as  that 
derived  by  Guderley  (reference  40)  by  an  entirely  different  (hodo- 
graph)  method  for  a  somewhat  differently  shaped  10-percent-thick 
airfoil.  Second,  the  value  of  c^  ^  =.0.017  given  at  MQ  =  1  by 

the  lower  transonic  formula  (equation  (103))  was  slightly  different 
from  the  value  0.011  obtained  at  M0  =  1  from  the  upper  transonic 

formula  (104).  This  slight  difference  is  due  to  r  approaching  a 
different  value,  2,  from  s( — *1)  as  M0  — >1.  Inasmuch  as  this 

difference  in  c^  ^  is  not  great,  and  in  view  of  the  considerable 

differences  in  procedure  on  the  two  sides  of  sonic  speed,  the  small 
difference  in  ^  is  to  be  regarded  not  as  a  defect  but  rather 

as  indicating  the  consistency  of  the  present  approach.  The  occur¬ 
rence  of  a  drag  maximum  at  M0  =  1  rather  than  somewhere  in  the 

vicinity  is  probably  due  to  the  approximations  involved  in  the 
present  calculation,  especially  that  which  assumes  that  the  terminal, 
shock  starts  at  the  trailing  edge. 

As  the  free- stream  Mach  number  increases  from  sonic  speed,  the 
drag  coefficient  decreases  because  of  the  decreasing  over-all  pres¬ 
sure  level.  As.  the  shock  attaches  and  the  field  of  flow  becomes 
completely  supersonic,  the  drag  finally  approaches  the  Ackeret  rate 

of  variation  with  Mach  number,  1//\|mo2-1.  Several  drag  coefficients 

have  been  calculated  in  the  completely  supersonic  range  by  using  the 
exact  Prandtl-Meyer  pressure  distribution  with  initial  conditions  at 
the  leading  edge  given  by  the  exact  oblique- shock  relations  of  refer¬ 
ence  49.  The  difference  between  the  first  such  point,  c^  =  0.051 

at  M0  =  1.50,  at  which  sonic  speed  occurs  behind  the  attached  shock 
at  the  leading  edge,  and  the  point  c^  =  0.065  at  M0  =  1.31  by  the 
present  method,  indicates  again  the  approximation  introduced  into  the 
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present  results  primarily  by  the  transonic  approximation  jl  <<  1 
(that  is,  setting  all  the  T  values  equal  to  1.2 )  and,  probably 
to  a  lesser  extent,  by  the  small-perturbation  approximation 

A  <<  l. 


CONCLUDING  REMARKS 

The  small-perturbation  form  of  an  integral  method  has  been  used, 
with  the  aid  of  simplifying  assumptions  indicated  by  experiment  or 
other  theory,  to  calculate  approximate  flow  patterns  for  a  symmetric 
biconvex  airfoil  through  the  transonic  speed  range.  The  wholly  sub¬ 
sonic  symmetric  flow  field  develops  symmetric  continuous  locally 
supersonic  regions  when  the  subsonic  free- stream  Mach  number  exceeds 
a  definite  (lower  critical)  value.  Above  an  equally  definite  (poten¬ 
tial  limit)  free-stream  Mach  number,  no  continuous  symmetric  solu¬ 
tion  exists  by  the  present  method.  Above  the  potential-limit  Mach 
number,  the  method  yields  both  a  physically  unreal  symmetric  solution 
containing  forbidden  regions  bounded  by  limiting  lines  and  a  physi¬ 
cally  real  asymmetric  solution  containing  an  asymmetric  locally  super¬ 
sonic  region.  Termination  of  this  locally  supersonic  region  on  the 
downstream  side  by  a  shock  required  additional  hypotheses  concerning 
the  shock.  These  hypotheses  were  quite  unrelated  to  the  limiting- 
line  phenomenon.  It  was  found  possible  to  join  the  asymmetric  solu¬ 
tion  continuously  to  the  asymmetric  solution  that  exists  for  slightly 
supersonic  free-stream  Mach  numbers.  The  main  results  of  the  method, 
such  as  the  variation  of  pressure  drag  with  Mach  number,  relative 
constancy  of  local  Mach  number  ahead  of  the  terminal  shock  in  the 
transonic  range,  and  so  forth,  agree  qualitatively  with  experimental 
data.  More  extensive  experimental  data  are  needed  in  the  form  of 
pressure  distributions,  wake  surveys,  and  terminal- shock  surveys  on 
airfoils  of  simple  shape  under  isolated  airfoil  conditions  at  high 
Reynolds  numbers.  The  present  method  is  probably  simple  and  flex¬ 
ible  enough  to  enable  suitable  alterations  and  extensions  to  be  made, 
whether  for  more  rigorous  solutions  of  the  cases  treated,  or  to 
more  general  situations,  such  as  airfoils  with  lift,  bodies  of 
revolution,  and  so  forth. 


National  Advisory  Committee  for  Aeronautics, 
Lewis  Flight  Propulsion  Laboratory, 

Cleveland,  Ohio,  December  19,  1949. 
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APPENDIX  A 


SMALL-PERTURBATION  APPROXIMATION  OF 
CONTINUITY  INTEGRAL 

The  exact  form  of  the  equation  of  continuity,  equation  (1), 
can  "be  written  as  follows  (see  fig.  14) : 


r?, 


nyo+6 


dY< 


pv  cos  0  4y 


(Al) 


where  5  is  a  possible  lateral  displacement  of  a  streamline  far 
from  the  airfoil.  Subtracting  y0-Y  from  both  sides  of  equa¬ 
tion  (Al),  as  in  the  section  COMPRESSIBLE  POTENTIAL  FLOW,  and  pro¬ 
ceeding  to  the  limit  yQ — >  oo  yields 

Y-S  b  (pv  cos  0-l)dy  (A2) 

JY 

An  expression  for  dy  in  terms  of  v  and  9  is  next  derived 
starting  from  the  differential  equations  of  motion  in  intrinsic 
form.  These  equations  are  (reference  35) 


Bv 

Bn 


-Cv 


(A3) 


where 


(A4) 


(A5) 


and  the  subscript  s  denotes  values  taken  along  a  streamline. 
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Consider  the  variation  of  velocity  v  and  stream  direction  0 
along  a  line  making  an  angle  aQ  vith  the  free-stream  direction 
(fig.  14).  If  z  denotes  distance  measured  along  this  line,  then, 
using  equations  (A3)  to  (A5), 

5v  /5v\  /50\  5s  5v  bn 
bz  =  \50/s  \dsjg  bz  +  bn  bz 

E!l.\  C  cos  (ao-0)  -  Cv  sin  (ao-0) 

bejg 


=  -Cv  sin  (ao-0) 


v(le^  COt  ^a°"0^  +  1 


bO  _  /50\  5s  bd  5n 
5z  \5a/s  bz  5n  bz 


rn  -C  cos  (ao-0)  -  C  ($L)  sin  (ao-0) 

v  I  be/ 


=  -C  sin  (ao-0) 


cot 


(A6) 


(A7) 


Equations  (A6)  and  (A7)  yield,  respectively,  for  variations  along  z, 

dv 


dz  sin  (ao-0) 


v  + 


5v 

59 


•j  cot  (aQ-0) 


(AS) 


dz  sin  (ao-0) 


d0 


cot  (ao-e)  +  <M  -11  (*) 
v  ^50/j 


(A9) 


If  the  angle  a0  is  now  taken  as  jt/2,  the  line  element 
dz  "becomes  dy,  and  equations  (A8)  and  (A9)  become 
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dy  cos  0  = 


dv 


C 


v  +  tan  0 


(A10) 


dy  cos  0  - 


C 


d0 


tan  0  + 


(M2-l)  /dv\ 

V  Ws 


(All) 


Substituting  equation  (A10),  the  definition  v  =  1+ v,  the  small- 
perturbation  expressions  cos  0=1-  02/ 2,  tan  0=0,  and  the 
expansion  equation  (10)  of  p  in  powers  of  v  into  equa¬ 
tion  (A2)  yields 

poo 

Y-5  = 

Jy 


llU-Tu2  -  £_-pV  qv 


0s 


dU 

C 


(A12) 


in  which 


(A13) 

(A14) 


The  next  higher  order  terms  in  the  bracket  of  equation  (A12)  would 
be  of  order  v^,  V02,  \iv^9vq3  .  The  order  of  magnitude  of  Qa 
is  v.  if  equation  (A12)  is  consistent  with  the  transonic  similarity 
rules,  then  02  ~  pv2.  The  two  terms  involving  0  in  equa¬ 
tion  (A12)  therefore  can  at  most  only  affect  the  function  rc(M0) 
but  not  its  limiting  value  rc(l)  =  T  =  (y+l)/2.  Thus  the  small- 
perturbation  equation  (4)  of  the  continuity  integral  is  correct  to 
the  order  v  in  the  subsonic  range  p  ~  1  and  to  the  order  V2 
with  rc  =  T  in  the  transonic  range  p  «1.  The  small-perturbation 
form  of  the  irrotationality  condition  (equation  (5))  is  seen  from 
equation  (A10)  to  be  correct  to  the  first  order  in  v  and  0. 


In  the  important  special  case  where  the  airfoil  has  fore-and- 
aft  symmetry,  equations  (4)  and  (5)  hold  exactly  at  the  midchord 
location,  for  here,  by  symmetry,  0=0  in  equations  (Al)  and  (A10). 
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In  the  neighborhood  of  such  a  location  and,  more  generally.  In  the 
neighborhood  of  any  0=0  location  on  an  airfoil,  it  may  therefore 
be  expected  that  the  essential  small-perturbation  approximations 
are  those  involving  the  velocity  increment  v,  as  in  equation  (10) 
rather  than  those  involving  the  flow  direction  0.  This  expectation 
is  the  justification  for  the  use  of  Tc  rather  than  r  in  the 
continuous -potential -flow  calculations  of  this  paper. 
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APPENDIX  B 


SMALL-PERTURBATION  APPROXIMATION  OP 


MOMENTUM  INTEGRAL 

The  momentum  principle'  in  integral  form  states  that  the  vector 
rate  of  flow  of  momentum  out  of  a  fixed  volume  of  space  equals  the 
vector-integrated  pressure  on  the  hounding  surface  of  the  volume . 
Applying  this  principle  to  the  two-dimensional  region  ABDE  (fig.  14) 
gives,  for  components  in  the  x -direction. 


-y0 

^yo+5 

^0+8  P 

(Po+Povo2)Ay  + 

JO  \- 

PoAy  = 

Yo  J 

(p+pv2  cos2  0)dy  + 

Y  Jo 

pady 


ED 


DB 


AB 


EOA 


(Bl) 

in  which  the  paths  of  integration  are  indicated  under  the  respective 
integral  signs.  Rearranging  the  terms  of  equation  (Bl)  as  with 
the  continuity  equation  ( Al) ,  proceeding  to  the  limit  y0 — >  co } 


denoting  by  p,  p,  and  v  fractions  < 
stream  quantities, and  noting  that  P0v0‘ 


>f  the  respective  free- 
'/p0  =  7M02  results  in 


Y-5 


n7, 


sJO 


7Mq2 


dY  + 


psr  +  (pv2  cos  2  6-1) 

bMo2 


Ay  (B2) 


Subtraction  of  the  continuity  equation  (A2)  then  yields  for 
the  pressure-drag  coefficient  c^  f  on  the  airfoil  from  the  leading 
edge  to  the  point  A  (including  a  factor  2  for  both  surfaces  of  the 
symmetric  airfoil). 


r»Y 


'd,f 


4 
7«o 


(p_-l)dy  *  -4 


Uo 


OY 


Eli—  +  pv  cos  0(v  cos  0-1)1  dy 

7Mq2  J 

(B3) 
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Substituting  the  definition  v  *  1+p,  the  expansion 
cos  0=1-  02/2  +  .  .  .,  the  expansions  for  p  and  p  as  power 
series  in  v  obtained  from  equation  (6),  and  the  irrotationality 
expression  (A10)  into  equation  (B3)  yields 


cd,f 


-  Up2 

2 


+  !r 


+  T  f  I  v2qvob 


dv 

~ 


(B4) 


where 


r 


m  = 


3 

2 


(B5) 


and  Uqs  is  defined  in  equation  (A14).  Equation  (B4)  shows  why 
the  momentum  integral  equation  (B2)(or  rather  the  derivative  of 
equation  (B2)  with  respect  to  Y)  is  inconvenient  to  use  along 
with  equation  (4)  in  solving  the  velocity-distribution  problem. 

The  flow-direction  term  02/ 2  in  equation  (B4)  is  of  the  same 
order  of  magnitude  as  the  velocity  term  -(|i/2)p2  and  so  cannot 
be  neglected  in  the  small-perturbation  case,  as  was  the  same  term 
in  the  continuity  integral  equation  (A12).  The  other  0  term  in 
equation  (B4)  is  of  order  and,  as  in  the  continuity  integral, 

at  worst  merely  modifies  the  function  rm(Mq)  but  not  its  limit¬ 
ing  value  at  sonic  speed  rm(l)  =  T  *  (l+r)/2.  Thus,  in  this  paper, 
of  the  three  first  integrals  of  the  differential  equations  of 
motion  -  that  is,  conservation  of  energy,  mass,  and  momentum  -  the 
first  two  (the  first  as  Bernoulli's  equation)  were  used  in  calcu¬ 
lating  the  basic  flow  pattern,  and  the  x-component  of  the  third 
integral,  on  the  basis  of  the  calculated  flow  pattern,  was  used 
to  determine  a  portion  of  the  pressure  drag.  The  lift  and  the 
moment  of  the  forces  acting  on  different  portions  of  the  airfoil 
could  similarly  be  calculated  from  conditions  in  the  flow  field 
by  the  y-component  momentum  integral  and  the  angular  momentum 
integral,  respectively. 


For  the  calculation  of  the  pressure  drag  on  the  forward 
portion  of  the  biconvex  airfoil  in  the  lower  transonic  range,  the 
forward  portion  extended  to  midchord.  Hence  the  two  terms  involving 
0  in  equation  (B4)  were  considered  as  vanishing  to  a  sufficient 
approximation.  The  resulting  integration  of  equation  (B4)  using 
the  curvature  function  (equation  (33))  then  yields  equation  (103). 

In  the  upper  transonic  range,  the  forward  portion  of  the  airfoil 
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extended  to  the  0  =  6V  location  and  6  was  held  constant  at 
this  value  in  the  integration  of  equation  (B4).  The  term  involving 
v@3  in  equation  (B4)  vanishes  in  the  limiting  case  Mq  »  1  or 
rm  =  T,  for  which  the  calculations  in  this  range  were  made.  Inte¬ 
gration  of  equation  (B4),  using  the  curvature  function  and  other 
relations  pertaining  to  the  upper  transonic  range  (see  section 
UPPER  TRANSONIC  PLOW  WITH  DETACHED  SHOCK)  leads  to  equation  (104). 
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APPENDIX  C 

TRANSONIC  SMALL-PERTURBATION  APPROXIMATION  OF 
OBLIQUE-SHOCK  RELATIONS 

An  element  of  oblique-shock  front,  inclined  at  an  angle  a 
to  a  free -stream  flow  at  Mach  number  Mo  (fig.  9)  produces  a  flow 
deflection  9  given  by  (reference  49) 


(T±\ 

4 

\7+y 

«  +  JL  -  COS  dQ 

L(7-1)M02  J 

sin  2 a  tan  (o-9) 


(01) 


The  corresponding  velocity  ratio  v  Just  downstream  of  the  element 
is 


v 


cos  g 
cos  ( o~9) 


(C2) 


Elimination  of  o  between  equations  (Cl)  and  (C2)  yields  the 
equation  of  the  shock  polar.  A  relation  between  the  slope 
(dv/daJsb  of  the  tangent  to  the  shock  polar,  that  is,  the  rate 
of  change  of  velocity  ratio  with  respect  to  stream  direction  along 
the  downstream  side  of  the  shock  and  the  similar  rate  of  change 
along  a  streamline  (dv/de)3  at  the  same  downstream  point  of  the 
shock,  is  obtained  by  dividing  equation  (B8)  by  equation  (A9)  and 
replacing  aQ  by  <j  ;  thus. 


4 


(b,\  T  *  (til cot  (0  -9) 

'  coMc-*,+^(|)g 


(C3) 


A  similar  relation  for  the  rate  of  change  of  velocity  vith  stream 
direction  in  the  y -direct  ion  at  a  point  on  the  downstream  side  of 
the  shock  element  is  given  by  setting  a  *  x/2  in  (C3), 


v  +  tan  9 


(C4) 
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The  ratio  of  the  local  curvature  Csjj  =  -  (do/dz)^  of  a  shock 
front  to  the  curvature  of  the  streamline  just  downstream  of  the 
shock  element  is  obtained  by  dividing  equation  (A8)  by  (do)^, 
setting  a0  =  a,  and  going  to  the  transonic  small-perturbation 
limit.  Thus, 


_  p 

csh  r(-ti  +  COS  a) 

—  *  - 1 - £ -  <C14> 

cos  o  (3|a-7  cos  a) 

The  conditions  at  various  points  on  the  shock  polar  in  the  small - 
perturbation  transonic  range  can  now  be  determined  by  the  preceding 
equations.  Crocco's  point,  for  example,  is  defined  by  (bv/bO)3  coo. 
Equations  (C12)  and  (C14)  show  that  this  condition  implies  zero 
streamline  curvature  if  the  shock  front  itself  has  finite  curvature, 
which  is  usually  the  case  for  a  detached  shock  wave.  Crocco's  con¬ 
dition  yields,  by  equation  (C12), 

cos2  o  «  3/7  JT  (C15) 

for  the  inclination  of  the  shock  front  in  the  physical  plane  at 
Crocco's  point.  Using  equation  (C15),  the  various  flow  quantities 
V,  9 ,  and  the  local  Mach  number  M  just  downstream  of  the  shock 
element  can  be  determined  from  equations  (C6),  (C7),  and  (27) 

with  rM  *=  r  • 


The  derivative  (be/bv)y  (equation  (C13))  is  infinite  at 
Crocco's  point,  zero  at  the  sonic  point  (defined  by  M  *  1),  and 
zero  also  at  a  subsonic  point  of  the  shock  polar  defined  by 
cos  a  =>  l/5  |i.  At  this  last  point,  the  downstream  streamline 
starts  out  at  constant  pressure,  inasmuch  as  here  (dc/de)s  =  0  by 
equation  (C12). 

The  properties  of  various  points  on  the  shock  polar  in  the 
small-perturbation  transonic  limit,  determined  in  the  manner  indi¬ 
cated,  are  listed  in  similarity  form  in  table  II.  Equations  (64) 
and  (68)  to  (70)  are  obtained  from  the  values  of  Py  /jT  given  in 
the  table.  The  velocity  increment  v  in  the  table  is  (vg/v^)  -  1, 
and  hence  equals  A^-Ag  in  the  section  LOWER  TRANSONIC  ASYMMETRIC 
FLOW  WITH  TERMINAL  SHOCK.  Similarly,  {T  in  table  I  is  i_i+2rAx 
in  the  same  section. 

Further  discussion  of  the  shock  polar  in  the  transonic  small- 
perturbation  limit  is  given  by  Guderley  (reference  38),  who  shows, 
in  particular,  that  from  Crocco's  point  to  the  point  of  maximum 
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*  deflection  on  the  shock  polar,  the  oblique  shock  in  the  physical 

plane  is  attached  to  the  vertex  of  a  wedge  with  curved  sides  and 
has  infinite  curvature  at  the  vertex  (Csh/C  =<»);  reference  38 
also  shows  that  the  point  of  maximum  deflection  on  the  shock  polar 
corresponds  to  the  point  of  shock  detachment  from  the  vertex, 
which  had  "been  doubted  by  Crocco  (reference  46)  and  more  recently 
by  Thomas  (reference  50). 


TABU 

5  II 

Defini¬ 

tion 

l/&0\ 

2 

:  cos^a 

cos  a 

Tv 

re 

P/2 

r1/3!) 

M2-l 

U 

\T 

e2/3 

Mach 

wave 

v  *  0 

-1 

i 

1 

0 

0 

0 

1 

Sonic 

point 

EB 

0 

1/2 

0.707 

-0.500 

0.354 

-1 

0 

Crocco’ s 
point 

(b0\  rs 

(  T““  )  ass  0 

Ws 

CO 

3/7 

0.655 

-0.571 

0.374 

-0.143 

Maximum 

deflec¬ 

tion 

(b e\ 

Wsh  =  0 

-i/a/F 

1/3 

0.577 

-0.667 

0.385 

-1.260 

-0.300 

Constant 

pressure 

fbv\ 

hrr  *  0 

\50/  g 

0 

1/5 

0.447 

-0.800 

0.358 

-1.587 

-0.600 

Normal 

shock 

6  m  0 

CO 

0 

0 

m 

0 

OO 

-1 

APPENDIX  D 


TERMINAL-SHOCK  CONDITION  IN  LOWER 
TRANSONIC  RANGE 

The  condition  that  for  finite  perturbation  of  the  flow  the 
terminal  shock  in  the  lower  transonic  range  he  such  that  its  back 
pressure  is  as  close  as  it  can  get  to  free-stream  pressure  may  be 
regarded  as  a  rough  statement,  for  a  specific  purpose,  of  the 
theorem  that  the  equations  of  motion  can  be  derived  from  the  vari¬ 
ational  principle 


(Po-P)  dx  dy  =  0  (Dl) 

where  p  is  local  pressure  at  a  point  (x,y)  and  pG  is  the  free- 
stream  pressure  at  infinity.  This  theorem  has  been  established  for 
nonviscous  flow  fields  containing  curved  shocks  (references  16,  43, 
and  51  and  unpublished  results  of  Chi -Teh  Wang) .  It  is  assumed  here 
to  hold  at  high  Reynolds  number  in  the  lower  transonic  range  More¬ 
over,  the  extremum  value  of  the  integral  in  equation  (Dl)  is  assumed 
to  be  the  same  as  for  irrotational  flow,  namely,  a  minimum. 

Consider  a  stream  tube  of  flow  Just  outside  the  boundary 
layer  (fig.  15).  In  this  stream  tube,  the  flow  accelerates  isen- 
tropically  from  free-stream  pressure  and  free-stream  Mach  number 
at  station  0  to  a  supersonic  condition  at  station  1,  reverts 
through  a  normal  shock  to  subsonic  conditions  at  station  2,  and 
returns  ultimately  to  free-stream  pressure,  but  not  to  free-stream 
velocity,  at  station  3.  The  flow  from  station  2  to  station  3  need 
not  be  assumed  lsentroplc.  If  the  normal  shock  is  assumed  to  occur 
in  a  diverging  portion  of  the  stream  tube,  as  is  Indicated  by 
stability  considerations  (reference  52),  then  the  maximum  velocity, 
or  minimum  pressure,  on  the  airfoil  occurs  Just  ahead  of  the  shock 
at  station  1.  If  the  shock  does  not  form  in  a  diverging  portion 
of  the  stream  tube,  then  the  pressure  coefficient  to  be  derived  is 
that  Just  ahead  of  the  shock,  not  the  minimum  in  the  stream  tube. 
The  small -perturbation  assumption  is  not  made. 

The  shock  back  pressure  .  P2  is  now  expressed  as  a  fraction 
of  free-stream  pressure  po.  By  the  isentropic  relation,  the  ratio 
of  the  pressure  Just  ahead  of  the  shock  p^  to  the  free-stream 
pressure  pD  is  given  by 
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The  Ideal  pressure  ratio  across  the  shock  is  given  by  the  Rankine- 
Hugoniot  relation 


27  ,,  2  7-1 

7+1  Ml  "  7+1 


(D3) 


As  in  the  section  on  lower  transonic  asymmetric  flow,  to  allow  for 
the  possibility  that  the  pressure  rise  across  the  actual  shock 
region  of  finite  extent  is  different,  for  whatever  reason,  from 
the  pressure  rise  across  an  ideal  normal  shock,  a  shock-pressure - 
drop  coefficient  is  introduced:  (The  desirability  of  considering 
the  influence  of  the  shock-boundary -layer  interaction  in  this  con¬ 
nection  was  pointed  out  to  the  author  by  A.  Silverstein. ) 


Aps  Aps 

sc  .  „  as  . 

1^2  w  2 

2^1t1  2PlMl 


(D4) 


such  that  the  actual  pressure  ratio  across  the  shock  is  given  by 


Equations  (D2)  to  (D5)  combine  to  give  for  the  actual  shock  back 
pressure  ratio 


For  given  M0  and  cs,  this  expression  is  found  to  have  a  maxi¬ 
mum  at  a  value  of  given  by 
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M 


7+3  -  ca(  r-i-l) 


2  -  c8(7+1)/2 

The  maximum  shock  hack  pressure  ratio  is  accordingly 


(D7) 


7 


7 


7-1 

2 


7-1 

|  (7+1-Cg ) 


4-cs(7+l) 
(7+1) (7+1-Cg) 


7-1 


(D8) 


It  is  now  inferred,  from  the  condition  that  shock  hack  pressure  he 
as  close  as  possible  to  free-stream  pressure,  that  if  (P2/P0)mfa:  <1> 

then  the  shock  Mach  number  Mq  is  given  by  equation  (D7),  whereas 
if  (p2/pQ)mflT>  1,  the  shock  Mach  number  Mq  is  given  by  equa¬ 
tion  (D6)  with  P2/P0  =  !• 

The  ideal-shock  case  cs  =  0  is  discussed  first.  Here  the 
free-stream  Mach  number  Mq,  at  which  (P2/p0)majc  *  1#  Is  0.78. 
Hence,  by  equation  (D7),  the  shock  Mach  number  Mq  for  free-stream 
Mach  numbers  Mq  <  0.78  is  given  by 


Mi.  1.483  (7*  1.4)  (D9) 

For  free-stream  Mach  numbers  in  the  range  0.78  <  Mq  <1,  equa¬ 
tion  (D6)  with  P2/P0  =  1  an(l  c8  =  0  gives  a  relation  between 
Mq  and  Mo.  The  preceding  relations  between  Mq  and  Mq  can  be 
converted  into  relations  between  the  shock  pressure  coefficient 


P  s  Pl"Po  2  ^P1 

1  |povo2  7M02  \p0  / 


(DIO) 


and  Mq  by  equation  (D2).  The  result  is  shown  in  figure  16. 

Shown  for  comparison  are  the  pressure-coefficient  curves  for  sonic 
velocity  (Mq  =  1  in  equation  (D2))  and  for  perfect  vacuum  (Mq  =  » 
in  equation  (D2)). 
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Included  in  figure  16  are  some  experimental  data  given  in 
references  44  and  53.  The  data  of  reference  53  consist  of  the 
minimum  pressures  measured  on  the  wings  of  various  airplanes  in 
both  low-speed  supercritical  flight  (stalled  pull-ups)  and  in  high¬ 
speed  supercritical  flight  (dive  pull-outs).  Of  the  six  circle 
data  points  in  figure  16,  the  three  points  at  free-stream  Mach 
numbers  of  0.42,  0.445,  and  0.645  represent  minimums  of  all  the 
pressures  measured  on  the  wings  of  the  P-47C-1,  the  XP-51,  and  the 
SB2C-1  airplanes,  respectively,  in  stalled  pull-ups.  The  cluster 
of  three  points  at  approximately  M0  =  0.73  represents  minimum 
pressures  selected  in  reference  53  from  the  pressure  distributions 
measured  on  the  wings  of  the  XP-51  airplane  in  successive  dive 
pull-outs.  A  further  dive  pull-out  point  for  the  XF2A-2  airplane 
coincides  with  the  middle  one  of  this  cluster  of  three  points. 

The  data  of  reference  44  were  obtained  on  a  NACA  4412  airfoil 
spanning  the  Langley  24 -inch  high-speed  closed  circular  tunnel. 

The  variation  of  minimum  pressure  on  the  upper  surface  with  Mach 
number  at  an  angle  of  attack  of  the  airfoil  of  l052.5t  was  taken 
as  representative  of  their  results.  Further  wind-tunnel  data  of 
the  same  type  are  given  in  reference  24. 

Comparison  of  the  experimental  data  with  the  theoretical  limit 
curve  shows  that  the  limit  pressures  indicated  by  the  theory  are, 
in  fact,  approached  or  reached  but  not  exceeded.  The  leveling  off 
of  the  pressure  distribution  of  reference  44  somewhat  below  a  local 
Mach  number  of  1.5  may  be  due  to  low  Reynolds  number,  wind-tunnel 
effect,  or  both.  Thus  the  concept  of  a  limit  curve  derived  from  a 
shock  condition  independent  of  Reynolds  number  and  airfoil  shape 
appears  to  have  some  significance,  at  least  in  the  limit  of  large 
Reynolds  numbers  (single  normal  terminal  shock),  isolated  airfoils, 
and  for  an  airfoil  shape  that  permits  high  local  Mach  number  super¬ 
sonic  regions  to  develop  at  low  free-stream  Mach  numbers. 

Consider  now  the  effect  of  shock  pressure-drop  coefficient  cs. 

As  noted  in  the  section  LOWER  TRANSONIC  ASYMMETRIC  FLOW  WITH 
TERMINAL  SHOCK,  values  of  c3  in  the  approximate  range  0.1  to  0.5 
occur  experimentally.  Assume  for  definiteness  a  value  c-  =  0.33. 

The  effect  of  this  value  of  c3  in  equations  (D6)  to  (D8)  is  to 
increase  the  shock  Mach  number  Mp  yielding  maximum  P2/P0  to  1.5 
and  to  shift  the  free-stream  Mach  number  M0  at  which  (P2/p0)mAT  «  1 

to  M0  =  1.  Hence,  the  portion  ABCD  of  the  cs  =  0  limit  curve 
(fig.  16)  would  be  eliminated  and  replaced  by  the  curve  BE,  slightly 
shifted  to  correspond  to  Mi  =  1.5  instead  of  1.483.  The  intro¬ 
duction  of  a  shock  pressure-drop  coefficient  cs  thus  yields  prac¬ 
tically  the  same  value  of  limit  shock  Mach  number,  and  hence  mini¬ 
mum  pressure,  as  when  c8  =  0  in  the  range  Mq  <0.78.  It  may  be 
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noted  that  the  value  c3  =  0.33  yields  the  pressure  ratio 
Pg/p0  =  0.527  in  equation  (D6),  that  is,  the  pressure  ratio  cor¬ 
responding  to  locally  sonic  velocity  in  isentropic  flow,  which 
provides  an  interpretation  for  the  statement  occasionally  made 
that  sonic  velocity  occurs  behind  the  terminal  shock. 

The  shock  condition  as  thus  far  developed  applies  to  finite 
and,  in  fact,  quite  large  velocity  perturbations  and  where  1-M0 
can,  and  possibly  should,  be  of  the  order  1.  [The  transition  to 
the  transonic  small-perturbation  case  is  equivalent  to  replacing 
the  branch  AB  of  the  limit  curve  (fig.  16)  (or  its  modification 
for  O3  /  0)  by  the  tangent  to  AB  at  the  point  A.  This  fact  is 
seen  from  equation  (67),  which  is  the  limiting  form  of  equa¬ 
tion  (D6)  for  |i  =  1-Mq^ — >0  and  M^-l — >0.  For  constant  ca 
near  Mq  =  1  and  by  the  shock  back  pressure  condition  P2  =  0  in 
this  range,  equation  (67)  yields  -2A.-^(=  P-jJ  as  a  linear  function 
of  1-Mq,  tangent  to  equation  (D6)  at  Mq  «*  1. 

As  remarked  in  the  section  LOWER  TRANSONIC  ASYMMETRIC  FLOW  WITH 
TERMINAL  SHOCK,  application  of  the  shock  condition  Pg/Po  =  1  or 
Pg  =  0  in  equation  (67)  implies  approximately  zero  pressure  change 
along  the  airfoil  downstream  of  the  shock.  If  this  condition 
exists  regardless  of  airfoil  shape,  a  severe  separation  of  the 
boundary  layer  must  occur  in  traversing  the  shock.  Similarly,  for 
finite  perturbations  at  lower  free-stream  Mach  numbers,  the  present 
analysis  applies  only  if  the  boundary  layer  can  so  vary  in  thick¬ 
ness  as  to  provide  sufficient  flexibility  in  the  airfoil  boundary 
condition  for  the  terminal  shock  to  be  free  to  stabilize  at  dif¬ 
ferent  local  Mach  numbers  along  the  airfoil.  Only  then  does  the 
possibility  arise  of  choosing  one  of  these  local  Mach  numbers  in 
accordance  with  equation  (Dl)  for  shock  location,  independently  of 
airfoil  shape. 

With  such  freedom  of  shock  location,  equation  (Dl)  also  Indi¬ 
cates  that  the  terminal  shock  will  form  in  a  diverging  portion  of 
stream  tube,  as  this  provides  less  variation  of  static  pressure 
than  if  the  shock  formed  (at  the  same  local  Mach  number)  in  a  con¬ 
verging  portion  of  the  stream  tube,  which  would  necessarily  be  pre¬ 
ceded  by  a  diverging  part.  Where  substantial  pressure  change  does 
occur  downstream  of  the  shock,  the  terminal- shock  condition  given 
herein  requires  modification;  for  example,  in  the  small-perturbation 
case  the  P2  term  in  equation  (67)  would  have  to  be  retained. 

Finally,  it  is  of  interest  to  note  that  th9  shock  Mach  number 
given  by  equation  (D9)  Is  also  that  for  which  the  streamline  curva¬ 
ture  vanishes  just  downstream  of  a  normal -shock  element  of  finite 
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curvature  (references  18  and '43).  This  property  is  derived  by 
finite -perturbation  considerations  analogous  to  those  of  appendix  C 
thus  no  close  connection  with  the  present  derivation  of  equa¬ 
tion  (D9)  appears  to  exist.  Inasmuch  as  a  separating  boundary 
layer  traversing  a  shock  should,  however,  tend  to  produce  zero 
streamline  curvature  downstream  of  the  shock  element,  this  property 
of  equation  (D9)  may  also  help  to  explain  terminal-shock  stabili¬ 
zation  at  approximately  Mq  =  1.5. 


NACA  TN  2130 


APPENDIX  E 

SYMBOLS 


The  more  important  symbols  used  in  this  report  are  as  follows; 
streamline  curvature. 


-  W 

\5s Jt 


chord  length,  taken  as  unity 

D(drag)/|-  p0v02c,  section  drag  coefficient 
1  2 

Aps/—  p-jv^  ,  shock-pressure  -drop  coefficient 

U 


drag 

distance  of  normal  element  of  detached  shock  from  leading  edge 

-2Ki  =  (1-Mo2)/(tD2/3,  transonic  similarity  parameter, 
(equations  (12c)  and  (84)) 

I^A/t2/3,  velocity  parameter,  (equations  (58b)  and  (85)) 

local  Mach  number,  ratio  of  local  velocity  to  local  velocity 
of  sound 

(p-p0)/|p0v02,  pressure  coefficient 

pressure  (fraction  of  free-stream  value) 

curvature  function  parameters,  (equations  (33)  and  (74)) 

1  +  A,  resultant  velocity  at  airfoil  (fraction  of  free-stream 
value) 

1  +  D,  resultant  velocity  (fraction  of  free-stream  value) 

i  (1+E),  chordvise  location,  distance  in  free-stream  direction 


ordinate  of  airfoil 

lateral  distance  perpendicular  to  free-stream  direction 
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z  =  u/A,  (distance  in  appendix  C) 
2 


r  =  1  + 


n  = 
p 
0 
T 

X  = 


-g-  -  Mq2>  (equation  (lie)) 


r.  =  (7+l)/2,  (equation  (lid)) 


+  zirlio*)’ 


r  jyj  =  Mo  [  1  +  Mq  j  ,  (equation  (28a) ) 

y  *  ratio  of  specific  heats,  1.4 
8  streamline  displacement  at  infinity,  figure  2 


angle  "between  local  stream  direction  and  free-stream 
direction,  figure  2 


\  (equation  (37a)) 

\3-r/  n 


-  U  =  1-Ho 

density  (fraction  of  free-stream  value) 
angle  "between  shock  element  and  free-stream  direction 
thickness  ratio,  ratio  of  maximum  thickness  to  chord  length 

1 —  YC&,  (equation  (37b)) 

(3-r)2  u3 

~  of  the  order  of  magnitude  of 

Subscripts : 

a  at  airfoil 

e  calculated  sonic  point  of  Prandtl -Meyer  distribution 

i  incompressible  or  low  speed 

I  potential  limit 
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o  free-stream  condition 

w  sonic  point  of  oblique  shock 

1,2  upstream  and  downstream  of  terminal  shock,  respectively 
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Figure  1.  -  Schematic  superposition  of  three  transonic  flow 
patterns  about  symmetric  airfoil* 
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Figure  2#  -  Symmetric  compressible  potential  flow# 
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Figure  4.  -  Velocity  distributions  onT=0.1  Kaplan  section  in  continuous  potential  flow, 
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Figure  5#  -  Effect  of  Mach  number  on  maximum  velocity  increment  of 
various  T  =  0.1  sections#  x  a  0*5# 
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(b)  Symmetric  biconvex  section. 


Figure  5 


#  -  Continued.  Effect  of  Mach  number  on  maximum 
increment  of  various  T  s  0.1  sections.  x  =  0.5. 
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(c)  Elliptic  section# 


Figure  5#  -  Concluded,  Effect  of  Mach  number  on  maximum  velocity 
increment  of  various  T  =  0,1  sections,  x  s  0.5, 


Jacobian  =  pv°  3(v,e)/3(  4>#¥) 
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Figure  6.  -  Jacobian  of  transformation  from  physical  to  hodograph 
plane.  Symmetric  biconvex  airfoil,  T=  0.1. 
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Figure  7.  -  Velocity  distributions  on  T  =0.1  symmetric  biconvex 
airfoil  in  subsonic  and  transonic  range# 


Terminal- shock  location,  x^ 
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Free-stream  Mach  number,  M0 


Figure  8  •  —  Terminal- shock  location  on  T=0.1  symmetric 
biconvex  airfoil  in  lower  transonic  range.  Limiting 
calculation  for  T— >0,  Mo  — *1*  Plotted  by  transonic 
similarity  rules  for  T  =  0.1, 
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Figure  9.  -  Upper  transonic  or  detached  shock  range. 
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Figure  10.  -  Detached  ^shock  location  relative  to  i  =  0.1  symmetric 
biconvex  airfoil  in  upper  transonic  range*  Limiting  calculation 
for  T— >0,  M0— >1*  Plotted  by  transonic  similarity  rules  for 

r  =  o*i. 
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(a)  x  s  0.5. 

Figure  12.  -  Local  Mach  number  on  symmetric  biconvex  airfoil  In 
transonic  range.  Limiting  calculation  for  T — »0,  Mo — *>1. 
Plotted  by  transonic  similarity  rules  for  t  s  0.1. 
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(b)  z  ■  0.8* 


Figure  12*  -  Concluded.  Local  Mach  number  on  symmetric  biconvex  air¬ 
foil  in  transonic  range.  Limiting  calculation  for  M),  Uq — >1. 
Plotted  by  transonic  similarity  miles  for  'f  s  0.1. 
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